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(57) Abstract 

An NMR data set (310, 31 1) which is acquired at (218) with applied phase encoding and read-out gradients at (136) is cor- 
rected to reduce motion artifacts and increase image sharpness. The acquired NMR data set (310, 31 1) is examined at (100) to de- 
tect bulk displacements of the object being imaged and phase displacements caused by motion. Tnis information is employed to 
produce correction operators (330, 331, 337) which are applied to the NMR image data set. One or more navigator signals (340) 
may also be acquired at (218) during the scan to produce an NMR data set from which the corrective operators can more readily 
be derived. 
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REDUCING MOTION ARTIFACTS IN NMR IMAGES 

« 

Cross Reference 

This application is a continuation-in-part of U.S. 
Application Serial No. 07/276,168, filed November 23, 
1988, 

5 Background of the Invention 

The field of the invention is nuclear magnetic 
resonance imaging methods and systems . More 
particularly, the invention relates to a method for 
reducing image artifacts caused by flow and motion. 
10 Any nucleus which possesses a magnetic moment 

attempts to align itself with the direction of the 
magnetic field in which it is located. In doing so, 
however, the nucleus precesses around this direction at 
a characteristic angular frequency (Larmor frequency) 
15 which is dependent on the strength of the magnetic field 
and on the properties of the specific nuclear species 
(the magnetogyric constant y of the nucleus) . Nuclei 
which exhibit this phenomena are referred to herein as 
"spins". 

20 When a substance such as human tissue is subjected 

to a uniform magnetic field (polarizing field B z ) , the 
individual magnetic moments of the spins in the tissue 
attempt to align with this polarizing field, but precess 
about it in random order at their characteristic Larmor 
25 frequency. A net magnetic moment M 2 is produced in the 



direction of the polarizing field, but the randomly 
oriented magnetic components in the perpendi collar , or 
transverse, plane (x-y plane) cancel one another* If, 
however, the substance, or tissue, is subjected to a 
magnetic field (excitation field Bx) which is in the x-y 
plane and which is near the Larmor frequency, the net 
aligned moment, M 2r may be rotated, or "tipped", into 
the z-y plane to produce a net transverse magnetic 
moment Mx, which is rotating, or spinning, in the x-y 
plane at the Larmor frequency. The degree to which the 
net magnetic moment M z is tipped, and hence the 
magnitude of the net transverse magnetic moment Mx 
depends primarily on the length of time and the 
magnitude of the applied excitation field Bx. 

The practical value of this phenomenon resides in 
the signal which is emitted by the excited spins after 
the excitation signal Bx is terminated. In simple 
systems the excited spin induce an oscillating sine wave 
signal in a receiving coil. The frequency of this 
signal is the Larmor frequency, and its initial 
amplitude, Ao, is determined by the magnitude of the 
transverse magnetic moment Mi. The amplitude,, A, of the. 
emission signal decays in an. exponential fashion with 
time, t: 

A - A 0 e- t / T * 2 

The decay constant 1/T*2 depends on the homogeneity, 
of the magnetic field and on Ta, which is referred to as 
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the "spin^spin relaxation" constant, or the "transverse 
relaxation" constant* The T2 constant is inversely 
proportional to the exponential rate at which the 
aligned precession of the spins would dephase after 
5 removal of the excitation signal Bi in a perfectly 
homogeneous field* 

Another important factor which contributes to the 
amplitude A of the NMR signal is referred to as the 
spin-lattice relaxation process which is characterized 

10 by the time constant Ti . This is also called the 

longitudinal relaxation process as it describes the 
recovery of the net magnetic moment M to its equilibrium 
value along the axis of magnetic polarization (z) . The 
Ti time constant is longer than T2, much longer in most 

15 substances of medical interest. 

The NMR. measurements of particular relevance to the 
present invention are called "pulsed NMR measurements". 
Such NMR measurements are divided into a period of 
excitation and a period of signal emission. Such 

20 measurements are performed in a cyclic manner in which 

the NMR measurement is repeated many times to accumulate 
different data during each cycle or to make the same 
measurement at different locations in the subject. A 
wide variety of preparative excitation techniques are 

25 known which involve the application of one or more 
excitation pulses <Bi) of varying magnitude and 
duration. Such excitation pulses may have a narrow 
frequency spectrum (selective excitation pulse) , or they 
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may have a broad frequency spectrum (nonselective 
excitation pulse) which produces transverse 
magnetization Mi over a range of resonant frequencies. 
The prior art is replete with excitation techniques that 
. 5 are designed to take advantage of particular NMR 

phenomena and which overcome particular problems in the 
NMR measurement process. The present invention may be 
used with any of these pulse sequences* 

When utilizing NMR to produce images/ a technique 
10 is employed to obtain NMR signals from specific 

locations in the subject. Typically, the region which 
is to be imaged (region of interest) is scanned by a 
sequence of NMR measurement cycles which vary according 
to the particular localization method being used. The 
15 resulting set of received NMR signals are digitized and 
processed to reconstruct the image using one of many 
. well known reconstruction techniques - To perform such a 
scan, it is, of course, necessary to elicit NMR signals 
from specific locations in the subject . This is 
20 accomplished by employing magnetic fields (G X r Gy/ and 
G 2 ) which have the same direction as the polarizing 
field. Bo, but which have a gradient along the respective 
x, y and z axes. By controlling the strength of these 
gradients during each NMR cycle, the spatial 
25 distribution of spin excitation can be controlled and 
the location of the resulting NMR signals can be 
identified. 



WO 90/05920 



PCT/US89/04945 



5 

NMR data for constructing images can be collected 
using one of many available techniques, such as multiple 
angle projection reconstruction and Fourier transform 
(FT). Typically, such techniques comprise a pulse 
5 sequence made up of a plurality of sequentially 

implemented views. Each view may include one or more 
NMR experiments, each of which comprises at least an RF 
excitation pulse and a magnetic field gradient pulse to 
encode spatial information into the resulting NMR 
10 signal. As is well known, the NMR signal may be a free 
indication decay (FID) or, preferably, a spin-echo 
signal. 

The preferred embodiments of the invention will be 
described in detail with reference to a variant of the 
15 well known FT technique, which is frequently referred to 
as "spin-warp". The spin-warp technique is discussed in 
an article entitled "Spin Warp NMR Imaging and 
Applications to Human Whole-Body Imaging" by W.A. 
Edelstein et al., Physic ill MediciPfi and Biolocry. Vol. 
20 25, pp. 751-756 (1980) . 

Briefly, the spin-warp technique employs a variable 
amplitude phase encoding magnetic field gradient pulse 
prior to the acquisition of NMR spin-echo signals to 
phase encode spatial information in the direction of 
25 this gradient. In a two-dimensional implementation 

(2DFT) , for example, spatial information is encoded in 
one direction by applying a phase encoding gradient (G y ) 
along that direction, and then a spin-echo signal is 
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acquired iri the presence of a read-out magnetic field 
gradient (G x ) in a direction orthogonal to the phase 
encoding direction. .The read-out gradient present 
during- the spin-echo acquisition encodes spatial 
information in the orthogonal direction- In a typical 
2DFT pulse sequence, the magnitude of the phase encoding 
gradient pulse G y is incremented (AG y ) in the sequence 
of views that are acquired during the scan to produce a 
set of NMR data from which an entire image can be 
reconstructed. 

Object motion during the acquisition of NMR image 
data produces both blurring and "ghosts" in the phase- 
encoded direction* Ghosts are particularly apparent 
when the motion is periodic, or nearly so. For most 
physiological motion each view of the NMR signal is 
acquired in a period short -enough that the object may be 
considered stationary during the acquisition window. In 
such case the blurring and ghosting is due to the 
inconsistent appearance of the object from view to view. 
Motion that changes the appearance between views such as 
that produced by a patient moving , by the respiration or 
the cardiac cycle/ or by peristalsis, is referred to 
hereinafter as "view-to-view motion". Motion may also 
change the. amplitude and phase of the NMR signal as it 
evolves during the pulse sequence and such motion is 
referred to hereinafter as "in-view motion". 

Both blurring and ghosting can be reduced if the 
data acquisition is synchronized with the functional 
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cycle of the object to reduce view-to-view motion. This 
method is known as gated NMR scanning, and its objective 
is to acquire NMR data at the same point during 
successive functional cycles so that the object "looks" 
5 the same in each view. The drawback of gating is that 
NMR data may be acquired only during a small fraction of 
the object's functional cycle, and even when the 
shortest acceptable pulse sequence is employed, the 
gating technique can significantly lengthen the data 

10 acquisition time. Some of these methods are disclosed 
in U.S. Patent Nos. 4,751,462; 4,567,893 and 4,663,591. 
None of them have proven entirely satisfactory because 
they either depend upon perfectly periodic motion, or 
they increase the scan time significantly, or they 

15 produce low signal-to-noise images. 

Several NMR pulse sequences have been proposed to 
either desensitize the NMR measurement to the phase 
perturbations caused by flowing spins as described in 
U.S. Patent No. 4,728,890, or to sensitize it to flow in 

20 such a manner that the effects of flow can properly be 
separated from the reconstructed images as described in 
U.S. Patent No. RE 32,701. None- of these methods have 
proven entirely satisfactory, either from a performance 
standpoint, or because of their adverse impact on scan 

25 time or the type of NMR measurements that may be 
performed. 

In our prior U,S. Patent No. 4,715,383, we disclose 
a method for reducing motion and. flow artifacts in NMR 
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images. While this method substantially improves NMR 
images by suppressing artifacts caused by spins outside 
the region of interest/ it does not correct for the 
motion artifacts produced by spins located, inside the 
5 region of interest . 

All of the prior methods for reducing motion and 
flow artifacts focus on the data acquisition procedure. 
They change the NMR pulse sequence itself , they change 
the order in which the pulse sequences in a scan are 

10 executed, . or they synchronize the execution of the pulse 
sequence with the motion of the subject under study. 
Their objective is to produce a set of NMR data which is 
minimally affected by flow and motion and which can, 
therefore, be used to . construct a clear, ghost-free 

15 image. 

. Summary of the TnrantH on 

The present invention is a method for reducing 
motion and flow artifacts in an NMR image by correcting 
the set of NMR data which has been acquired during a 

20 scan to remove the effects of motion and flow before the 
image is reconstructed. More specifically, the 
invention includes transforming . the NMR data set to 
create a hybrid^space data array; producing a correction 
data array using- the data in the hybrid-space data 

25 array;, and applying the data in the correction data 

array to the NMR- data set produced by the NMR system to 
reduce flow, and motion artifacts in the image which is 
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reconstructed from the NMR data set. The correction 
data array which is produced in accordance with one 
aspect of the present invention corrects for view-to- 
view motion artifacts and the correction data array 
5 which is produced in accordance with another aspect of 
the invention corrects for in-view motion and flow 
artifacts. 

A general object of the invention is to provide a 
motion and flow artifact correction method which can be 

10 employed after the NMR data has been acquired. The 
present invention enables the NMR data set to be 
corrected retrospective of its acquisition and, 
there fore f it may be used in addition to any motion or 
flow artifact suppression techniques which have been 

15 employed in the past. To the extent that the acquired 
NMR data is affected by motion or flow, the present 
invention will detect it and automatically correct the 
NMR data so that the affect is reduced or eliminated. 
Another aspect of the invention is to acquire 

20 "navigator" NMR data along with the usual image NMR data 
within the same pulse sequence. The navigator NMR data 
enables the corrections for view-to-view and in-view 
motion artifacts to be made more accurately. A 
navigator NMR signal is produced in each pulse sequence 

25 along with the image NMR signal and a data set is 

acquired for both. The corrective values are determined 
using the navigator data set and the corrections are 
made to the image data set. 
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A more specific object of the invention is to 
correct the NMR data set for view-to-view artifacts 
caused by motion or -f low irx any direction. The 
navigator signail can be acquired in the presence of a 
5 read-out magnetic field gradient which is oriented in 

any direction. The shift corrections which are produced 
according to the present invention, will correct for 
errors caused by view-to-view motion or flow in the 
direction of the navigator signal read-out gradient. By 

10 acquiring more than one navigator signal in the pulse 

sequence in the presence of read-out gradients oriented 
in respective different directions,, shift corrections 
are produced which correct for view-to-view motion in 
the corresponding directions. For example, shift 

15 corrections can be made for motion along both the x axis 
and the y axis of the NMR system. 

Yet another specific object of the invention is to 
correct, the NMR data set for in-view artifacts caused by 
motion orr flow in any direction. The phase corrections 

20 E which are produced according to the present invention 
offset artifact causing systematic noise produced by in- 
view motion or flow regardless of its direction. 

Still another aspect of the present invention is 
that both ih-view and view-to-view motion and flow 

25 artifact correction can be applied to the same NMR data 
set . Both the in- view correction data array and the 
view-to-view correction data array can be applied to the 
NMR data set. 
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The foregoing and other objects and advantages of 
the invention will appear from the following 
description. In the description/ reference is made to 
the accompanying drawings which form a part hereof, and 
5 in which there is shown by way of illustration a 

preferred embodiment of the invention. Such embodiment 
does not necessarily represent the full scope of the 
invention r however, and reference is made therefore to 
the claims herein for interpreting the scope of the 
10 invention. 

Brief Description of the Drawings 

Fig. 1 is a block diagram of an NMR system which 
employs the present invention; 

Fig. 2 is an electrical block diagram of the 
15 transceiver which forms part of the NMR system of Fig. 
i; 

Fig. 3 is a graphic representation of a 
conventional NMR pulse sequence used to acquire data to 
produce an image; 
20 Fig. 4 is a pictorial representation of how an 

image is reconstructed from NMR data acquired using the 
pulse sequence of Fig. 3; 

Fig. 5 is a pictorial representation of how 
correction values are calculated according to the 
25 present invention; 

Fig. 6 is a graphic representation of data in the 
modulus array of Fig. 5; 
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F iQ"- 7 is a graphic representation of the 
correlation process used to produce the shift correction 
array of Fig. 5; 

Figs.. 8a and 8b are graphic representations of data 
5 in the phase array of Fig. 5 and the process used to 
produce the phase correction array of Fig. 5; 

Fig. 9 is a graphic representation of an alternative 
pulse sequence used to acquire both image NMR data and 
navigator NMR data;: and 
10 Fig.. 10 is a plot of hybrid-space data which 

illustrates the effects of both stationary and moving 
spins on hybrid- space phase. 

DescrlPi- ion of the Ersfexrgd Embodiment 

Referring to Fig. l r there is shown in block 
diagram form the major components of a preferred NMR 

15 system which incorporates the present invention and 

which is sold by the General Electric Company under the 
trademark "SIGNA" . The overall operation of the system 
is under. the control of a host computer system generally 
designated 100 which includes a main computer 101 (a 

20 Data General .HP4QQ0)'* The computer 100 includes an 
interface 102 through which a plurality of computer 
peripheral devices, and other NMR system components are 
coupled to the main computer 101. Among the computer 
peripheral devices is a magnetic tape drive 104 which 

25 may be utilized under the direction of the main computer 
101 for archiving patient data and image data to tape. 
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Processed patient data may also be stored in an image 
disc storage device designated 110. An array processor 
106 is utilized for preprocessing acquired NMR data and 
for image reconstruction. The function of image 
processor 108 is to provide interactive image display 
manipulation such as magnification, image comparison, 
gray-scale adjustment and real time data display. The 
computer system 100 also includes a- means to store raw 
NMR data (i.e. before image construction) which employs 
a disc data storage system designated 112. An operator 
console 116 is also coupled to the main computer 101 by 
means of interface 102, and it provides the operator 
with the means to input data pertinent to a patient 
study as well as additional data necessary for proper 
NMR system operation, such as calibrating, initiating * 
and terminating scans. The operator console is also 
used to display images stored on disc or magnetic tape. 

The computer system 100 exercises control over the 
NMR system by means of a system control 118 and a 
gradient amplifier system 128. Under the direction of a 
stored program, the computer 100 communicates with 
system control 118 by means of a serial communication 
network 103 (such as the Ethernet network) in a manner 
well known to those skilled in the art. The system 
control 118 includes several subsystems such as a pulse 
control module (PCM) 120, a radio frequency transceiver 
122, a status control module (SCM) 124, and power 
supplies generally designated 126. The PCM 120 utilizes 



.14 

control signals generated under program control by main 
computer 101 to generate digital waveforms which control 
gradient coil excitation, as well as RF envelope 
waveforms utilized in the transceiver 122 for modulating 
the RF excitation pulses. The gradient waveforms are 
applied to the gradient amplifier system 128 which is 
comprised of G x , G y and. G z amplifiers 130, 132 and 134 , 
respectively. Each amplifier 130, 132 and 134 is 
utilized to excite a corresponding gradient coil in an 
assembly designated .136 which is part of a magnet 
assembly 146. When energized, . the gradient coils 
generate magnetic field gradients G x , G y and G 2 . 

The gradient magnetic fields are utilized in 
combination with radio frequency pulses generated by 
transceiver 122, RF amp 123 and RF coil 138 to encode 
spatial information into the NMR signals emanating from 
the region of the patient being studied. Waveforms and 
control signals provided by the pulse control module 120 
are utilised by the transceiver subsystem 122 for RF 
carrier modulation and mode control. In the transmit 
mode, the transmitter provides a radio frequency signal 
to an RF power amplifier 123 which then energizes RF 
coils 138 which are situated within main magnet assembly 
-146. The NMR signals radiated by the excited spin in 
the patient are sensed by the same or a different RF 
coil than is used for transmitting. The signals are 
detected, amplified, demodulated, filtered, and 
digitized in the receiver- section of the transceiver 
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122. The processed signals are transmitted to the main 
computer 101 by means of a dedicated, unidirectional, 
high-speed digital link 105 which links interface i02 
and transceiver 122. 
5 The PCM 120 and SCM 124 are independent subsystems 

both of which communicate with main computer 101, 
peripheral systems, such as patient positioning system 
152, as well as to one another by means of serial 
communications link 103. The PCM 120 and SCM 124 are 
10 each comprised of a 16-bit microprocessor (such as Intel 
8086) for processing commands from the main computer 
101. The SCM 124 includes means for acquiring 
information regarding patient cradle position, and the 
position of the moveable patient alignment light fan 
15 beam (not shown) . This information is used by main 

computer 101 to modify image display and reconstruction 
parameters. The SCM 124 also initiates functions such 
as actuation of the patient transport and alignment 
systems . 

20 The gradient coil assembly 136 and the RF transmit 

and receiver coils 138, are mounted within the bore of 
the magnet utilized to produce the polarizing magnetic 
field. The magnet forms a part of the main magnet 
assembly which includes the patient alignment system 

25 148, a shim coil power supply 140, and a main magnet 

power supply 142. The main power supply 412 is utilized 
to bring the polarizing field produced by the magnet to 
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the proper; operating strength of 1.5 Tesla and is then 
disconnected* . 

To minimize interference from external sources, the 
NMR system components comprised of the main magnet 
assembly, the gradient coil assembly, and the RF 
transmit and receiver coils, as well as the patient- 
handling devices, are enclosed in an RF shielded room 
generally designated 144. The shielding is generally 
provided by a copper or aluminum screen network which 
encloses the entire room. The screen network serves to 
contain the. RF signals generated by the system, while 
shielding the system from RF signals generated outside 
the roonu- 

Ref erring particularly to Figs. 1 and 2, the 
15 transceiver 122 includes components which produce the RF 
excitation field Bi through power amplifier 123 at a 
coil 138A and components which receive the resulting NMR 
signal- induced -in a coil 138B. The base, or carrier, 
frequency of the RF excitation field is produced by a 
20 frequency synthesizer 200 which receives a set of 

digital signals through the communications link 103 from 
the main computer 101. These digital signals indicate 
the frequency which is to be produced at an output 201 
at a resolution of one Hertz. This commanded RF carrier 
25 is applied to a modulator 202 where it is frequency and 
amplitude modulated in response to signals received 
through line 203,* and the resulting RF excitation signal 
is turned on and off in response to a control signal 
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which is received from the PCM 120 through line 204. 
The magnitude of the RF excitation pulse output through 
line 205 is attenuated by a transmit attenuator circuit 
206 which receives a digital signal from the main 
5 computer 101 through communications link 103. The 

attenuated RF excitation pulses are applied to the power 
amplifier 123 that drives the RF transmitter coil 138A. 

Referring still to Figs. 1 and 2, the NMR signal 
produced by the excited spin in the subject is picked up 
10 by the receiver coil 138B and applied to the input of a 
receiver 207. The receiver 207 amplifies the NMR- signal 
and this is attenuated by an amount determined by a 
digital attenuation signal received from the main 
computer 101 through link 103. The receiver 207 is also 
15 turned on and off by a signal through line 208 from the 
PCM 120 such that the NMR signal is acquired only over 
the time intervals required by the particular 
acquisition being performed. 

The received NMR signal is demodulated by a 
20 quadrature detector 209 to produce two signals I and Q 
that are coupled through anti-aliasing filters 216 and 
217 to a pair of analog to digital converters indicated 
collectively at 218. The quadrature detector 209 also 
receives an RF reference signal from a second frequency 
25 synthesizer 210 and this is employed by the quadrature 

detector 20 9 to sense the amplitude of that component of 
the NMR signal which is in phase with the transmitter RF 
carrier (I signal) and the amplitude of that component 
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of the NMR signal.which is in quadrature therewith (Q 
signal) . 

The X and Q components of the received NMR signal 
are continuously sampled and digitized by the A/D 
converter 218 at a sample rate of 64 kHz throughout the 
acquisition period. A set of 256 digital numbers are 
acquired for each land Q component of the NMR signal, 
and these digital, numbers are conveyed to the main 
computer 101 through the serial link 105. 

The NMR system of Fig. 1 performs a series of pulse 
sequences to collect sufficient NMR data to reconstruct 
an image. One such pulse sequence is shown in Fig. 3. 
This sequence performs a slice selection by applying a 
90° selective. RF excitation pulse 300 in the presence of 
a z axis gradient pulse 301 and its associated rephasing 
pulse 302. After an interval TEi/2, a 180° selective RF 
excitation pulse 303 is applied in the presence of 
another z axis gradient pulse 304 to refocus the 
transverse magnetization at the time TEi and produce an 
20 echo NMR signal 305. 

To position encode the echo NMR signal 305, an x . 
axis read-out gradient pulse 306 is applied during the 
acquisition of the NMR signal 305. The read-out 
gradient frequency encodes the NMR signal 305 in the 
well known manner. In addition, the echo NMR signal 305 
is position encoded along the y axis by a phase encoding 
gradient pulse 307. The phase encoding gradient pulse 
307 has one strength during each echo pulse sequence and 
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associated NMR echo signal 305 r and it is typically 
incremented in steps through 256 discrete strengths (- 
128 to +128) during the entire scan. As a result, each 
of the 256 NMR echo signals 305 acquired during the scan 
is uniquely phase encoded. 

It is, of course, usual practice to repeat the 
pulse sequence for each phase encoding gradient value 
one or more times and to combine the acquired NMR 
signals in some manner to improve signal-to-noise and to 
offset irregularities in the magnetic fields. In the 
following discussion, it is assumed that such techniques 
may be used to acquire the NMR data set which is to be 
corrected. 

Referring particularly to Fig. 4, the acquired NMR 
data is stored in the data disk 112 (Fig. 1) in the form 
of two 256 X 256 element arrays 310 and 311. The array 
310 contains the in-phase magnitude values I and the 
array 311 contains the quadrature values Q. Together 
these arrays 310 and 311 form an NMR image data set 
which defines the acquired image in what is referred to 
in the art as "k-space". 

To convert this k-space NMR data set into data 
which defines the image in real space (i.e. Cartesian 
coordinates) , a two step Fourier transformation is 
performed on the I and Q arrays 310 and 311. The 
transformation is performed first in the read-out 
direction which is the horizontal rows of the arrays 310 
and 311 to produce two 256 X 256 element arrays 312 and 



WO 90/05920 



PCT/US89/04945 



20 

313. The array 312 contains the in-phase data and is 
labeled I 1 , while the array 313 contains the quadrature 
data and is labeled Q f . The X* and Q' arrays 312 and 
313 define the acquired image in what is referred to in 
the art. as "hybrid-space 1 '- This first transformation of 
the acquired NMR data set is expressed mathematically as 
follows^ 

Sxy(k-space) => s^ty (hybrid- space) 

- 255 

S *'y =25? I Sxy e-i2«xV256 (1> 

The second transformation is performed in the phase 
encoding direction which is the vertical columns of the 
arrays 312 and 313 to produce two 256 X 256 element 
arrays 314 and 315. The array 314 contains the 
transformed in-phase values ..and is labeled T n r while the 
array. 315 contains the quadrature values and is labeled 
Q n . This second transformation may be expressed 
mathematically as follows: 

Sx'y (hybrid-space) => S x « y t (real space) 
255 

Sx 'y f ■- 25? X s x t y e-i2*yy'/25e (2> 

"y=rO 

The arrays 314 and 315 are a data set which defines the 
acquired image in real space and the elements, thereof 
are used to calculate the intensity values in a 256 X 
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256 element image array 316 in accordance with the 
following expression: 

lm xy - V(I"xy> 2 + <Q" X y> 2 (3) 

The 256 X 256 elements of the image array 316 are mapped 
to the main operator console 116 (Fig. 1) for display on 
a CRT screen. 

The above described NMR system and pulse sequence 
is representative of the current state of the art. The 
diagnostic quality of the image which is obtained is 
determined by the extent to which the acquired NMR 
signals are degraded by superimposed thermal 
(statistical) noise and systematic (artifact) noise 
(which includes unsharpness) . The sources of thermal 
noise are well understood. The ratio of NMR signal 
intensity to thermal noise is determined by such factors 
as polarizing magnetic field strength, RF receiver coil 
configuration, the particular pulse sequence used, and 
the amount of signal averaging used. Systematic noise, 
mainly resulting from physiological motion, may degrade 
the NMR image far beyond the fundamental limit set by 
the thermal noise. Indeed, the diagnostic quality of 
many NMR images is limited far more by motion artifact 
and other forms of systematic noise than by intrinsic 
thermal noise. However, in contrast to thermal noise, 
there is no inherent limit on the extent to which 
systematic noise can be reduced. 
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It is the reduction of systematic noise and 
unsharpness which the present invention addresses. 
Whereas prior techniques address this problem with 
improvements to the NMR system hardware (for example, 
.cardiac and respiratory gating circuits) or with 
improvements to the pulse sequence (for example r motion 
desensitizing gradient pulses) , the present invention 
typically addresses the problem retrospectively. More 
specif ically, it is a discovery of the present invention 
that systematic noise can be easily detected by 
examining the NMR data set in hybrid-space. Systematic 
noise caused by view-to-view motion as well as 
systematic noise caused by in-view motion can be 
detected. Corrective operators can then be calculated 
and used to eliminate the systematic noise from the NMR 
image data set. The usual image reconstruction process 
can then be performed, using the corrected NMR image data 
set to produce an image which is substantially free of 
artifacts caused by motion. 

. Referring to Fig.. 5,. the first step is to produce 
the hybrid-space I 1 and Q* arrays 312 and 313 as 
described, above (equation 1) . The image data set itself 
can be used to produce the hybrid-space arrays 312 and 
313 , or as will .be described in more detail below, 
separate NMR data produced by a navigator NMR signal 
within the same pulse sequence can be used. In either 
case, corrective values for both view-to-view and in- 
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view systematic noise can be calculated from the hybrid- 
space arrays. 

Referring to Figs, 5 and 6, it has been discovered 
that view-to-view systematic noise can be detected in 
5 the magnitude information contained in the hybrid-space 
arrays 312 and 313. Accordingly, the next step in the 
process is to produce a 256 X 256 modulus array (M) 320. 
This is accomplished by calculating each element Mxy of 
the modulus array 320 from the corresponding elements 
10 I'xy and Q'xy of the hybrid-space arrays 312 and 313. 

Mxy = Vd'xy) 2 + (Q'xy) 2 < 4 > 

where: x = array column number, 1—256; and 
y = array row number, 1-256. 

It is a discovery of the present invention that 
15 view-to-view motion in the direction of the associated 
read-out gradient can be seen as a shifting left or 
right of the modulus data in each row of the modulus 
array 320. This is illustrated in Fig. 6, where each 
graph is the modulus values in a horizontal row of the 
20 modulus array 320 plotted as a function of its column 
number in the array 320. Although each graph is 
slightly different due to the differences in the amount 
of phase encoding, the graphs do have significant peaks 
which should occur at the same column number of the 
25 array 320. This is illustrated by the dashed lines 325 
and 32 6 through the peaks 327 and 328 in the graph of 
the first row (y = 1). The corresponding peaks 327' and 
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328 1 in the graph of the second row (y = 2) are shifted 
to the left of the dashed lines 325 and 32 6 r and the 
corresponding peaks 327." and . 328" in the last row of 
data (y = 25 6> are shifted to the right. Examination of 
the remaining 509 rows of the modulus array 320 would 
show similar shifts in varying amounts which are caused 
by view-to-view movement of the subject. 

The next step in the process is to determine the 
amount r S, by which each row of modulus data must be 
shifted to bring it into correlation with a reference 
row. This cross correlation may be performed in a 
number of ways. In the preferred embodiment a reference 
row (M r > in the modulus array 320 is selected and it is 
cross correlated with each of the other 255 rows M y . 
This is done by determining how much each row M y must be 
shifted left or right to maximize the sum of the product 
of its elements and the corresponding elements in the 
reference row M r «. 

In other words, the following expression is 
calculated and stored: 
256 

X < M Y*> " (5) 

x=l 

The data on row M Y is then shifted one position and 
the process is repeated. These calculations are 
performed for M y shifted between -32 and +32 positions, 
although these can.be extended if necessary for larger 
displacements. A plot of an exemplary correlation curve 
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which results from these sixty-four calculations is 
shown in Fig. 7. It can be seen that the peak in this 
exemplary curve occurs when the row of modulus data- M y 
is shifted to the left a few positions. Accordingly, a 
5 corrective value S y is determined by finding the number 
of shifts needed to produce the peak in the correlation 
curve. A shift correction value S y is calculated for 
each row (y — 1 to 256) of the mod u l u s a rr a y 3 20 an d i s 
stored in a 1 X 25 6 element shift correction array 330. 

10 There are many other ways to correlate each row in the 
modulus array 320. 

A very significant reduction in view-to-view motion 
in the x axis can be achieved by applying the correction 
values S directly to the modulus array 320 and 

15 reconstructing an image from this corrected NMR data. 

This is accomplished by shifting the data in each row of 
the modulus 320 by the amount indicated by the 
corresponding element of the shift correction array 330. 
Then, the values in the I* and Q' hybrid-space arrays 

20 312 and 313 are calculated using the corrected modulus 
values and the known phase angle for each element. An 
image can then be reconstructed in the usual fashion 
from the corrected I' and Q' hybrid-space arrays 312 and 
313 as described above. 

25 While the " above-described correction procedure 

provides significant improvement in image quality, 
further improvements can be made. When the shift 
correction S is made, a small amount of phase error is 
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introduced into the NMR- data set. This is due to the 
phase rollover which occurs in any NMR system as a 
result of -asymmetry in the echo signal in. its 
acquisition window.. The signal represented by the 
elements in each column of the I 1 and Q 1 hybrid arrays 
312 and 313 contain the same amount of rollover phase, 
and the amount of rollover phase changes linearly at the 
rollover rate Kr as the arrays are traversed from left 
to right through its columns (x » 1 to 256) . Thus, when 
data in a row is shifted left or right to make . the S 
correction discussed above , the data moves into a column 
with a different phase value than the column of its 
origin. This rollover phase can be corrected and is a 
function of the "amount which the data was shifted: 

Or = Kr * s. 

The rollback, correction values for each row in 
hybrid-space is calculated to produce a 1 X 256 element 
rollback correction array 331. Each element in the 
array 331. indicates the phase correction which must be 
made to all elements in its corresponding row of the I' 
and Q' hybrid- space arrays 312 and 313. How these phase 
corrections are actually made will be discussed in more 
detail below. 

The correction made thus far accounts for view-to- 
view motion along the read-out x axis. As will be 
discussed in more detail below, similar corrections can 
be made along "the phase encoding y axis by using a 
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special navigator signal in the same pulse sequence as 
the image signal and which is acquired in the presence 
of a y axis read-out gradient pulse. 

The above-described corrections do not account for 
5 phase errors in the NMR data due to in-view motion or 
flow. Such errors occur because the spins are moving 
during the pulse sequence. It is a further discovery of 
the present invention that these phase errors can be 
detected in the k-space or the hybrid-space NMR data 
10 set. 

Referring to Figs. 5 and 8, the procedure begins 
again with the I f and Q' hybrid-space arrays 312 and 
313. This data set is used to calculate the phase O of 
the acquired NMR signal at each of the 256 sample times 
15 and at each of the 256 views, A 256 by 256 element 

phase array 335 is produced, with each of its elements 
having a value which is calculated as follows : 

<& xy = TAN-1 (Q'xy/I'xy) 

It is a discovery of the present invention that the 
2 0 calculated phase <I> should have the same value along any 
column (y =1 to 256) of the phase array 335. This is 
shown graphically in Fig. 8a where the calculated phase 
O for three rows of data has been plotted as a function 
of column number (x = 1 to 256) . For any given NMR data 
2 5 set which is produced without phase encoding gradients 
(i.e. a navigator signal), these plots will be 
substantially the same if there are no. phase errors . To 



WO 90/05920 



PCI7US89/04945 



10 



28 

the extent that the calculated phase <£ differs in any 
column, that difference represents a phase error which 
should be corrected. This is graphically illustrated in 
Fig. 8b where two of the phase plots in Fig. 8a are 
superimposed on one another to reveal differences in 
their values over the central region. These differences 
are illustrated by the arrows Oi and 0 2 . 

To correct for the phase error r therefore, a 256 by 
256 element phase correction array 337 is produced. 
This is accomplished by establishing one row (y - 1) of 
phase data in the O array 335 as a reference and then 
finding the difference between the value of one of its 
elements and the same element in each of the other rows 
{ y 2* to. 256) of the O array 335. This is repeated for 
each element C x - 1 to 256) of the reference row (y - 
1) and. the calculated difference values (AO F ) are stored 
in the corresponding location in the phase correction 
array 337. Thus r each value in the phase correction 
array 337 indicates the amount by which the phase of 
each element of the NMR data set should be corrected. 
This corrects for in-view motion in any direction (x, y 
or z) . How. this correction is actually made will be 
described in more detail below. 

While separate • phase correction values AO P are 
25 calculated for each of the 256 by 256 elements of the 
array 337 in the preferred embodiment, it should be 
apparent that a -Mess rigorous approach can also be used. 
If the in-view motion, is known to occur over only a 
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small segment of the x axis field of view, then the 
calculation of corrective values AOp may be limited to 

that segment. Furthermore, if the in-view motion is 
uniform over that x axis segment; then perhaps a single 
5 value AOp will suffice as a correction over the entire 

segment. Since these corrections are applied 
retrospectively to the acquired NMR data set, it is 
contemplated that correction variables such as this will 
be under operator control and the radiologist can 

10 manipulate the corrections to obtain the image he needs 
in minimal time. 

When calculating the phase corrections which are to 
be made, one must consider whether both moving and 
stationary spins are contributing to the measured phase 

15 values. Where all the spins are moving, no further 

correction is necessary. However, when stationary spins 
are contributing substantial signal, then the phase 
correction values can be further refined. 

The effects of NMR signal components produced by 

20 both stationary and moving spins is illustrated in Fig. 
10. The points on this plot represent the I 1 and Q f 
values for the 256 data samples in a single column of 
the hybrid-space data sets 312 and 313. These points 
define a circle which is illustrated by dashed line 401 

25 having a center 402. If no motion were present, all of 
the points in this ring would cluster at one location. 
On the other hand, if all the spins are in motion, then 
the ring 401 would be formed, but its center 402 would 
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be at the origin (Q'«0, I r =0>. In the illustrated 
example, the signal components due to stationary spins 
are represented by the vector A which offsets the center 
402 from the origin- The signal component due to moving 
spins is represented by the vector B and its magnitude 
determines the size of the circle 401. The values (<£} 
in the phase array 335 discussed above indicate the 
phase angle of the combined stationary and moving spins. 
To provide a more accurate correction for in-view 
motion, therefore, the values in the phase array 335 
should be altered to indicate the phase <B« of only the 
iitoving spins . 

This alteration of the phase array values is 
performed on one column at a time. First, the 256 data 
points in the column of the hybrid-space arrays 312 and 
313 are applied to a curve fitting program which 
determines the center 402 of the circle 401. The 
coordinate values of the center 402 are then subtracted 
from the values in that column of the respective hybrid- 
space arrays 312 and 313. These altered values are then 
used to produce the altered phase array 335, which, in 
turn, is used to produce the more accurate phase 
correction array 337 as described above. 

While the phase error detection method described 
above works well when applied to a hybrid-space phase 
array which is derived from data acquired from navigator 
NMR signals, less rigorous phase error detection methods 
are also . possible using the NMR image data set itself. 
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For example , the I and Q arrays 310 and 311 can be used 
to calculate a 256 X 256 phase array in k-space. While 
the phase values in this array may not be useful at its 
outer boundaries, the phase information at the center 
5 column (x=128) where the peak of the NMR echo signal 305 
is sampled is usually unambiguous • One element in this 
column is selected as the reference and all the other 
values are compared with it to produce a 1 X 256 element 
phase correction array. These correction values AOp may 

10 be applied to the NMR data set as will now be described. 

All of the corrections calculated according to the 
present invention can be made to the NMR raw image data 
set I and Q in k-space. In k-space all of the 
corrections are implemented as a rotation of the phase 

15 of each NMR signal sample (Sxy ~ ^xy + jQxy) • The 

corrected NMR data set S'xy ™ay thus be calculated as 
follows : 

S'xy eiA* T Sxy (6) 

This rotation is implemented in the I and Q arrays 
20 310 and 311 in accordance with the following 
expressions: 

I xy (corrected) = IxyCOsAOT - QxySinA^T 
Qxy (corrected) - IxySinAOT + QxyCosAO? 

The phase change AO? is the arithmetic sum of the 
25 phase changes needed to make all of the above-described 
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corrections. Specifically, the total phase correction 
AO T is calculated as follows: 

AOt — -(A<£ R + K * A<£ ? + AO x + AO y ) (7) 

This expression includes the rollback correction 
5 AOr front the array 331 (Fig. 5) and the phase correction 
A<£p from the array 337 (Fig. 5) . The phase correction 
AOp is multiplied by a conversion factor K, however , to 
account for any difference in phase shifts between the 
image. NMR signal and a navigator NMR signal which may be 

10 used to determine the phase shift correction- If the 
image NMR data set is used to calculate the phase 
corrections r this conversion factor is "one". 
Otherwise^ the value , of K is measured by comparing the 
phase shift produced in image data (with phase encoding 

15 gradient applied* and the phase shift produced in the 

navigator data. In the alternative, the value of K can 
be calculated . 

The phase rotations A<5 X and A^ y are made to 

correct for view-to-view motion and. flow effects along 
20 the respective x and y axes. These phase corrections 
are calculated from the shift values S in the shift 
correction array 330 (Fig. 5> : 

A<£ x « S x * (x-(N x -l>/2> * 2tt/N x (9> 

where: S x =* shift correction from array 330 

25 x = sample number (i.e. 0 through 255> 

N x = total number of samples during 
read-out (256) . 
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AO y - S y * <y-(N y -l)/2) * 27C/N y (10) 

where : Sy = shift correction calculated from 

navigator signal acquired in 

presence of a y axis magnetic field 

5 gradient 

y = phase encoding view number (0 to 
255) 

Ny « total number of phase encoding 
views (256) . 

10 It should be apparent that the correction is 

different for each element of the NMR image data arrays 
I and Q because the total phase correction A<E>t is 

different for each x and y location in these arrays. 

As indicated above , the corrective methods of the 
15 present invention are applied to the NMR data set which 
has been acquired for the purpose of producing the 
desired image. While this same image data set can also 
be used to derive the corrective values which are to be 
applied to it, an alternative procedure is to produce a 

20 separate NMR data set during the same acquisition scan 
from which the corrective values can more accurately be 
derived. This separate NMR data set is produced by 
acquiring one or more "navigator" NMR signals during 
each pulse sequence. The main distinction of these 

25 navigator NMR signals is that they are not subject to 

the application of incremented phase encoding gradients . 
In some cases , a constant phase encode gradient may be 
used to elicit two-dimensional investigatory 
information . 
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A navigator pulse sequence is shown in Fig* 9. As 
with, the conventional pulse sequence of Fig, 3, the 
spins are excited by a 90° selective RF excitation pulse 
300 in the presence of slice select gradient pulses .301 
5 and .302, and at a time interval TEi/2 later, the spins 
are subjected to a 180° selective RF excitation pulse 
303 in the presence of slice select gradient pulse 304. 
No phase encoding gradients are applied, however, and at. 
interval TEi an NMR echo signal 340 is acquired in the 

10 presence of x axis read-out gradient pulse 306. A 

second 180° RF excitation pulse 341 is then produced in 
the presence of a slice select gradient pulse 342, and 
the y axis phase encoding pulse 307 is then applied. 
The resulting NMR echo signal 305 produced at TE2 is 

15 then acquired in the presence of read-out gradient pulse 
306 f to build the NMR image data set as described above. 

The NMR data set which is created by acquiring the 
NMR navigator signal 340 is subject to substantially the 
same motion effects as the NMR image data set which is 

20 acquired at the same time. However, because the 
navigator signals 340 are not phase encoded, the 
magnitude of the navigator NMR signals 340 are not 
significantly diminished at the extremities of the field 
of view. Corrective values can, therefore, be more 

25 readily calculated using the NMR navigator signal data 
set rather than the NMR image data set. The corrective 
values described above are thus determined by using the 
acquired NMR navigator signal data set, and these 



WO 90/05920 



PCT/US89/04945 



35 

corrective values are then applied to the NMR image data 
set. The corrected NMR image data set is then used to 
reconstruct an image. 

Referring again to Fig. 5, the NMR navigator signal 
5 data set is transformed first along the read-out 

gradient direction to produce the hybrid-space I 1 and Q f 
arrays 312 and 313. This data is then used as described 
above to produce a shift correction array 330, a 
rollback correction array 331 and a phase correction 
10 array 337. These corrective values are applied to the 
NMR image data set 310 and 311 (Fig. 4) as described 
above, and the corrected NMR image data set is processed 
in the usual manner to produce the corrected image array 
316 (Fig. 4) . 

15 The shift corrections S correct for view-to-view 

motion in the direction of the read-out gradient. 
Accordingly, the NMR navigator signal 340 in Fig. 9 is 
acquired in the presence of the read-out gradient pulse 
306 which is directed along the x axis and the shift 

20 corrections will compensate for motion along the x axis. 
Alternatively, the NMR navigator signal 340 can also be 
acquired in the presence of a read-out gradient pulse 
which is directed along the phase encoding, or y axis, 
as shown by the dashed lines 343. In yet another 

25 alternative, the pulse sequence can be modified to 
produce two NMR navigator signals, with one being 
acquired during an x axis read-out gradient pulse and 
the other being acquired in the presence of a y axis 



read-out gradient pulse. Either one of the navigator 
signal NMR data sets can then be used to produce the 
phase correction array 337 (Fig. 5> , while both 
navigator signal NMR data sets are used to produce 
corresponding shift correction arrays 330. A rollback 
correction array 331 is not always required for shift 
corrections made in the phase encoding direction, y. 
Both the x and y shift corrections can be made in 
hybrid-space, or they may be converted to phase angle 
corrections: which are made to the- NMR image data set in 
Jc-space as described above. 

The two-dimensional shift correction data can also 
be derived using "a single navigator NMR signal acquired 
during the application of a read-out gradient. As 
previously described, the magnitude of the navigator NMR 
signal yields the correction values S along the 
frequency encoding direction (x> . In addition, however, 
the phase of this navigator NMR signal can be used to 
calculate the correction values along the phase encoding 
direction (y ) as f ol lows : 




(11) 



where: A0 py is the measured phase shift, A y is 



the amplitude of the applied phase 



encoding gradient, and N y is the number of 



phase encoding increments (i.e. views). 
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The correction values for use in the above equation (10) 
for the y direction is, therefore, given by the 
following when N y is 256: 

_ A» pYl 256, j 
Y 27C( j-128.5) 

5 where: j is the view number 

While the navigator NMR signal is an echo signal in 
the preferred embodiment , this is not a necessary 
requirement. The navigator NMR signal may also be 
produced as a gradient recalled signal or as the result 
10 of three-dimensional data acquisition pulse sequence. 

The preferred embodiment of the above-described 
invention is implemented by a Fortran program which is 
executed by the main computer 101 for correcting in-view 
motion and view-to-view motion in the direction of the x 
15 axis read-out gradient. A listing of this program is 
provided in Appendix A. 

While the corrective values which are derived 
according to the present invention are employed 
retrospectively to improve image quality, they can also 
20 be employed to alter the scan. This is particularly 
true of the shift values S which correct for view-to- 
view motion. More specifically, as each pulse sequence 
is executed and a row of raw NMR data is acquired for 
the image arrays 310 and 311, the row of data may be 
25 transformed to hybrid-space (equation (1)) , the modulus 
values calculated (equation (4)), and a correlation made 
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with the previously acquired row of NMR data (equation 
C5>) . The resulting shift value S can then be used to 
alter the operating conditions of the NMR system before 
the next pulse sequence is executed. For example, where 
the shift value corrects for motion in the y axis 
direction, the phase of the RF carrier signal produced 
by the transmit frequency synthesizer 200 (Fig. 2) is 
changed by the amount -A«S y (equation (9)), Similarly, a 
shift correction for motion along the slice select 
direction (z axis) can. be made by changing the frequency 
of the RF carrier produced by the transmit frequency 
synthesizer 200. Shift corrections for motion along the 
read-out direction (x axis) can be made by changing the 
• location of the field of view along the x axis. 
15 The present invention addresses the problems 

currently associated with magnetic resonance 
angiography, in addition, this invention is necessary 
for quantitative tissue characterization using magnetic 
resonance. The invention is not limited to Fourier 
image reconstructions (i.e. 2DFT or 3DFT) and may be 
used with line scanning and other projection 
reconstructions. 
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APPENDIX A 



This is the main MARCO program. It gats input from tha user, such things 
as which file to work with. It than calls tha routina DO MARCO that 
does tha actual work of gatting valuas from tha raw data~file. Thasa 
valuas includa thosa for both x diraetion and intraviav motion cor- 
rection. Thasa valuas ara than writing to a taxt fila to be rtad by the 
polo program. 



BEGINNING OF MARCO. F77 SOURCE FILE 



! COPYRIGHT - MAYO FOUNDATION FOR MEDICAL EDUCATION AND RESEARCH 

IMPLICIT NONE 1988 

CHARACTER* 128 SOX7RCE FILE ! Nama of source fila 

CHARACTER* 12 B ASCII ?ATH ! Nana of fila for ascii output 

INTEGER* 4 I ERR ~ i Error; flag 

INTEGER* 4 I CHAN 1 i Channel numbar for sourca fila 

LOGICAL*2 CHOPPER, SHIM, HNEX, STRIP, D3FFT, GOT_IT 
LOGICAL*2 EXORCIST, HI LO SORT, CHOPLETSON, NOAP, NO_SORT 
INTEGER* 2 SKIP920 

INTEGER* 4 I BYTES 1 

INTEGER* 2 BASELINE SKIP 
INTEGER* 2 FRAMES PfR SLICE 
INTEGER* 2 DATA BlOCXS 

INTEGER* 2 VECTORC ! Count ar into VECTOR, below 

INTEGER* 2 NEX 
INTEGER* 2 NAVS 
INTEGER* 2 NECHOES 
INTEGER* 2 FRAME SIZE 

INTEGER* 2 THE ECSOl ! Echo from which to gat displacement 

C " ! VECOTR 

INTEGER*4 DELTAS (1024) 
REAL PHASES ( 1024 > 

C 
C 

INTEGER* 2 VIEW FROM 
INTEGER* 2 STARTO 

INTEGER* 4 III, I 

LOGICAL*2 DO DELTAS, DO_RQLL. DO_SHIFT 

INTEGER* 2 NA? ECHO, IMAGE_ECHO 

CHARACTER ANS7 ANSI, ANS2 , ANS3 

INTEGER* 4 GET I 

INTEGER* 2 NEX~IN RAW 

INTEGER* 2 START BELTA 

INTEGER* 2 STOP BELTA 

INTEGER*2 FRAME DELTA 

INTEGER*2 START~ROLL 

INTEGER*2 STOP ROLL 

INTEGER* 2 PHASE COL 

REAL ANGLE ROLL 

INTEGER*2 PHASE FRAME 
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WRITEC6.I42) 
142 FORMAT ( ' 1 r , //////) 

III - 0 

£2 T «£ T " -FALSE. 

^Gorfrl i ( :TR& G0T - IT) • AND - < IX1 •"?• "» 
III - lit + i * 

WRITE (6, 3 010). 
3010 FORMAT ( 1 +ENTER SOURCE FILE 



C 

C. . . 

c. . . 



READ(*. ' f*> ') SOURCE_FILE 
.00 I- I , 128 - 

EN IF^S0URCE,FItEfI ! I).EQ.'.<40>-) SOORCE_FILE (I: X) »"<0>" 

IF^liiifgifoT TUT^ 1 ^' °' ^^"VTESl, I ERR) 
IF «i?S5 R I EQ * 21 > "HEN 

ELSE XT (IERR .EQ> 20) THEM 

"cONTJHW™ ***** 172 J ^ 1 ZERO *»WH FILENAME 
ELSE 

PRfOT^'^?7) IN 0PE * IN * SOURCE FILE, CODE s f > IERR 
END IF ' • - 

END IF 
END DO 

XF (IIX .GE. 1Q J GO TO 9000 

"TOm2?SlKTK£f -H"™ OF ^ RAW DATA FILE HEADER TO DE- 
TERMINE CERTAIN THINGS ABOUT THE RAW DATA FILE. 

GO^TO^OOO* 0 * FR ° M SE *« NG COLLECTION TYPE, CODE: MERR 
END IF, 

IF NEX°IN P RAW ) 1*°" ( - N6T ' N0AP >) THEN . - 

END XJ ~ 

DO r » 1 r 1024 
DELTAS (I) - -42 



WO 90/05920 



PCT/US89/04945 



41 



PHASES (I) - -42 
END DO 

DO_DELTAS - .TRUE. 
FRAME__DELTA - -42 

DO WHILE ( (FRAME_OELTA . LT. 0) .OR. . 

# IE£?RL?S TA - GT - FRAMES PER S LI CE * NEX_IN RAW) ) 

WRITE(6,4020)FKAMES_PER SLICE*NEX IN RAW ~ ~ 

i220 M FORMAT ('+WHICH OF THE '7l3, * FRAMES DO YC 
FERENCE?— — — <GT> • , $) 

FRAME_DELTA • INT2 (GET I{)) 
EHD DO ~ 



YOU WANT TO USE FOR THE RE 



4010 



4011 



START_DELTA - -42 

DO WHILE ((START DELTA .LT. 0) .OR, 

# (START^DELTA . GT. FRAME SIZE) ) 
WRITE(6,4010) 

FORMAT ( '+ WHERE DO YOU WANT TO START TO GET SELECTIVE DELTAS? ~ 

— -<GT> • , $ ) 

START DELTA - INT2 (GET I(V) 
END DO * 
IF (START DELTA . EQ. 0) THEN 
START DELTA - 1 
STOP DELTA • FRAME SIZE 
ELSE ~" • 
STOP DELTA - -42 

DO WHILE ( (STOP DELTA .LT. 0) .OR. 

# (STOP_pELTA .GT. FRAME SIZE)) 
WRITE(6,4011) ~ 

FORMAT ( 1 + WHERE DO YOU WANT TO STOP GETTING SELECTIVE DELTAS? 

<GT>» ,$) 

STOP DELTA - INT 2 (GET !()) 
END DO~ ~ 
IF (STOP DELTA .LT. START DELTA) THEN 
I - ST5P_DELTA 
STOP DELTA - START DELTA 
START DELTA - I " 
END IF ~ 
END IF 

DO_ROLL - .TRUE. 
DO_SHIFT - .TRUE. 
PHASE_FRAME - -42 

DO WHILE ( (PHASE FRAME .LT. 0) .OR. 
\ (PHASE_FRAME .GT. FRAMES PER SLICE * NEX IN RAW) ) 

WRITE (6 , 4030) FRAMES PER SLICE*NEX IN~RAW ~ " 

FORMAT ( ' +WHXCH OF THE ',12,' FRAMES DO YOU WANT TO USE FOR THE PH 
AS E REFERENCE — <GT>» ,$) 

PHAS ENFRAME - INT2 (GET I() ) 
END DO ~ 



4030 



4003 



NAV ECHO - 0 
IMAGE_ECHO - 0 

DO WHILE ((NAV ECHO .LT. 1) .OR. (NAV ECHO .GT. NECHOES) ) 
WRITE ( 6 , 4 0 0 3TNECH0ES 

FORMAT ( 1 +OF THE Ml,* ECHOES , WHICH IS THE NAVAGATION ECHO?- 
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<GT>' , S) 

NAV_BCH0 - INT2 (GET_I()) 



END 05 



*l^«;iSJai&- M V 1 ?"'"- ^^ ECHO •«* NECHOES) ) 
~^-_..-E2^^<GT>i T $f f -' 11 -'' ECtt0ES ' raiCH " THE IMAGE ECHO? 



end^ E - ECH0 * IOT2 <«*_i(j) 

IF (CHOPLETSON). THEN 

NEX - 1 
ELSE 

NEX - NAVS * 2 
END IF 

IF (HNEX) THEN 

^NT^'UNABLE TO PROCESS KAIiF NEX DATA.' 

GO TO 9000 
END IF 

1F wtS2I' tiiJ^SSSSSL,' 011 - CCHOPLETSON). .OR. (NOAP) ) ) THEN 

^S£*'*NOT CHOPPER, NOAP r OR CHOPLETSON, TOABLE'to PROCESS. 
GO TO 9000 

END IF 

I - 1 

DO^WHILE^Cfl .M. 122) .AND. (SOURCE_FIUE (III)- . NE • CHAR(O))) 
END DO 

IF f.I . LT. 122 ) THEN 

IF ((SOURCE_FILEf 1-4:1-1) .EQ. 1 . RAW • ) .OR. 
ST ^|§TO«^FILE(I-4:I-i) .£Q. r . rav • ) ) THEN 

ELSE 

STARTO - I" 
END IF * 

DoT= P tSHT?4 OUR ff/ IIEf ls ™»-« // ' * CNC ' 

^JPS*-********* * CHAR(O) 
END DO 

ELSE 

£5 X 2 T ** t£RROR ' SOURCE FILENAME IS TOO LONG. 1 
GO TO 9000 
END IF 

IF^IERR^NE*© ) ^THEN^^^^~^™ ' STATUS " ' FRESH ' < IOSTAT-IERR) 

P* 1 ****' 'ERROR DETECTED DURING FILE OPEN/ CODE : MeRR 
GO TO 9000 
END IF 

IF (CHOPPER .AND. NOAP) THEN 
. WRlTE(6 r 3040) 

ELSE IF (CHOPPER .AND. .NOT. NOAP) THEN 
WRITE ( €j r 3 050 ) 
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END IF 

WRITE(6, 3060) FRAMES PER SLICE, FRAME SIZE 
WRITE (6, 3 070) NECHOES, NE3? 

3040 FORMAT (• SOURCE FILE: CHOPPER AND NOAP 1 ) 
3 050 FORMAT (' SOURCE FILE: CHOPPER AND AP* ) 

m5£2!.' PER SLICE • 1 , I3,T32, 'FRAME SIZE - ',13) 

3070 FORMAT ( ECHOES - ',I3,T32,'NEX - M3,/> 

CALL DO_MARCO (DO DELTAS, DO ROLL, DO SHIFT , FRAME DELTA, 

# START_DEETA,STOP_d£LTA,STAET ROLL, STOP~ROLL, ANGLE ROLL, 
2 PH^ COL, PH£SE FRAME , DELTASTPHASES , NEX_IN RAW , NECHOES 
J FRAME_SIZE,FRAMES_PER SLICE, DATA BLOCKS? ■ 

# ICHAN1 , NAV ECHO, IMAGE~ECHO, IERR) ~ 
VECTORC - FRAMES PER SLICE * NEX IN RAW 
IF (IERR .HE. 0) THEN ~ ~ 

PRINT* , "ERROR FROM DO MARCO, CODE- 1 , IERR 
GO TO 9000 
END IF 

WRITE (42^56J^pO — pELT^^FRAME^pELTA 2 START^DELTA. STOP_ DELTA , 

# 

356 FORMAT 
DO I - 

WRITE (42, 4242) I, DELTAS (I) , PHASES (I) 
4242 FORMAT(I3,I6,F20.15) 1 ' 

END DO 

CLOSE (42 , IOSTAT-IERR) 
9000 CONTINUE ! ERROR EXIT 
END 




**** gjjj) 0F MARCO. F7 7 SOURCE FILE ***** 
BEGINNING OF DO MARCO. F77 SOURCE FILE 



COPYRIGHT - MAVO FOUNDATION FOR MEDICAL EDUCATION AND RESEARCH - 1988 

SUBROUTINE DO_MARCO ( DO_DELTAS , DO_ROLL, DO SHIFT, FRAME DELTA , 

# START_DELTA, STOP DELTA, START ROLE,STOP ROLL, ANGLE ROLL, 

# PHASE_COL,PHASE_?RAME, DELTASTPHASES, NEX IN RAW , NECHOES , 

# FRAME_SIZE, FRAMES PER_S LICE, DATA B LOCKS 7l CHAN 1 , NAV ECHO, 

# IMAGE ECHO, IERR) ~ \ 
IMPLICIT NONE" 

INTEGER* 4 IERR J ERROR FLAG 

INTEGER* 4 DELTAS (*) ! ARRAY OF DELTA SHIFTS 

REAL PHASES (*) I ARRAY OF PHASE SHIFTS 
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c. 
c. 
c. 



c 
c 
c 



REAL FI9 
PARAMETER (FI9 



2-827435388) ! PI * 0.9 



OBTAINED • *«" BOOLEANS THAT DETERMINE WHAT VALUES ARE TO BE 

LO§I§U?i E-fiS?* 8 : ! T ° GET TOE DELTA SHIFTS 
lSIiS£*I 5n-?2^ T0 ^ THE ROLL BACK ANGLE - 

LOGICAL*2 DO^SHIFT I TO GET THE PHASE SHIFTS 

INTEGERS FRAME DELTA 
INTEGER* 2 START~DELTA 
INTEGER* 2 STOP EELTA 
INTEGER* 2 START ROLL 
INTEGER*2' STOP 50LL 
INTEGER* 2 PHASE COL 
INTEGER* 2 NEX IS RAW 
INTEGER* 2 PHASE FRAME 
REAL ANGLE ROLL"" 



INTEGER* 2 BUFFER(1S384) 



INTEGER*2 £BUF'(512)- 
INTEGERS STEP 
INTEGER** LN 
INTEGER** I, J, K, L2 r 
INTEGER*4 ICHAN1 
INTEGER** GET I 
REAL ROLLERl4,512) t ^ 
REAL Exiyi, Exi, Eyi r Exi2 

INTEGER* 2 FRAMES PER SLICE 
INTEGER* 2 DATA BlOCKS 
INTEGER* 2 ST0P~BLOCK 
INTEGER* 2 FRAME SIZE 
INTEGER* 2 NECHOES 
INTEGER* 2 NAVS 
INTEGER* 2 FOFF 

INTEGER*2 FRAME FROM 
INTEGER* 2 NEX FROM 
INTEGER*2 CHOPPER FAC 
INTEGER* 2 FHASE_FRAHE FROM 
INTEGER* 2 PHASE_NEX^F5oM 

INTEGER* 2 LOW CON 
PARAMETER (LOW^CON ■ -32) 
INTEGER* 2 HIGH CON 
PARAMETER (HIGH~C0N 
INTEGER* 2 ONE ~ 
PARAMETER (ONE - I) 
INTEGER* 2 TWO 
PARAMETER ( TWO - 2) 

REAL" MOD MAX 
REAL MOD225 



£°5 SEEING IN VIEWS, IF THERE 
^^512 SAMPLES IN A FRAME, FOUR NEX, 
ECHOES, 16384 WORDS WOULD BE 

HEAD 



L, M 



REAL, IMAG, MOD, ANGLE 



32) 
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REAL MOD_10 
INTEGER* 4 MAXOV1 
INTEGER* 4 B C 
INTEGER* 2 VECTORC 
INTEGER* 2 PHASEC 
REAL RR 

rIaL VIEW1 (-256:256) 

REAL VIEWN(-256:256) 

REAL THE_MODS (-256:256) 

REAL MOD MEANS (LOW_CON : HIGH_CON) 

INTEGER*? LOW HALF 

INTEGER* 2 HIGH_HALF 

INTEGER* 2 NAV ECHO 

INTEGER* 2 IMA5E_ECHO 

REAL PHASE_REF 

COMPLEX CI (512) 
REAL Bl(2,512) 
EQUIVALENCE (CI , Bl) 

DO I « 1 , FRAME_SIZE 

ROLLER(l,D - 0.0 

ROLLER(2,I) - 0.0 
END DO 

LOW HALF - 1-FRAME SIZE/2 
HIGH_ HALF » FRAME_SIZE/2 

IF (NEX IN RAW . EQ . 1 ) JHE N 
FRAME~FR5M - FRAME.DELTA-l 

PHAS E^FRAME^ FROM - PHASE_FRAME-1 

phase~nex_f5om ■ i 

END IF 

WRITE (6 , 2020) 

IF (DO_DELTAS) THEN 
K - 1 

DO WHILE (K • LT*. START_DELTA) 

ZBUF(K) - 0 

K - K + 1 
END DO 

DO WHILE (K .LE. STOPJDELTA) 

ZBUF(K) - X 

K - K + 1 
END DO 

DO WHILE (K .LE. FRAME_SIZ£) 
ZBUF(K) - 0 
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K m K + 1 

END DO 

DO I - -256 , 25S 

VIEWl(I) - 0.0 

VIEWN(I) - O.O 
END DO 

END IF 1 DO_DELTAS 

VECTORC « 0 
PKASEC ■» Q 

FOFP - FRAME_SIZE * 2 

IF (FRAMEJSIZE .EQ. 512) THEN 

STEP » * * NECHOES * NEX IN RAW 

LN Sfr " "™ — 

ELSE IF (FRAME SIZE .EQ. 256) THEN 

STEP m 2 * NECHOES * NEX IN RAW 

LN m 8 . - ~ ~ 

ELSE 

I ERR « 4201 

RETURN 
END IF 

STOP^BLOCK « 28 + STEP * FRAMES_PER_SLICE 

XF (STEP * FRAMES_PER_SLICE .GT. DATA BLOCKS) THEN 

PRINT* * ERROR , NOT ENOUGH DATA BLOCKS. ' 

I ERR - 420a 

RETURN 
END IF - 

IF (DO_ROLL) THEN 
DO' I « 1 0 , -1 

CALL LIBSRDB (28+STEP* (FRAME SIZE/ 2— I) , STEP, BUFFER, I CHAN 1, 1 ERR) 
IF(IERR.NE.O) THEN ~ 

IF ( (IERR .NE. -1) .AND. XIERR .NE. 24)) THEN 
PRINT*, »•-■■-. 

PRINT *, 'ERROR DETECTED WHEN READING BLOCK NUMBER r * ,28+(S 
TEP*FRAKE_SIZE/2-I) r * CODE: ' , IERR 
ELSE 

PRINT* , ' ■ 

PRINT* , 1 ERROR , END OF FILE. 1 

mg. i f 

RETURN 
END IF 

DO L2 m o r NEX IN_RAW - 1 

B_C - FOFF * T(IKAGE_ECHO-I) f L2 * NECHOES) 

IF (KOD4(L2 r 2) .EQ. I) THEN 

CHOPPER FAC » -1 
ELSE - ~ 

. CHOPPER FAC « 1 
END IF ~ 
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L - FRAME SIZE / 2 + 1 

DO J - 1+B C , FRAME SIZE+B C , 2 

C1(L) « CMFLX( BUFFER (J )"* CHOPPER_FAC , 
BUFFER (J+l) * CHOPPER FAC) 

L « L + 1 
END DO 

L - 1 

DO J - FRAME SIZE+l+B C , FRAME SIZE*2+B_C , 2 
C1(L) - CM?LX( BUFFER ( J ) * CHOPPER_FAC , 
BUFFER (J+l) * CHOPPER FAC) 

L - L + 1 
END DO 

CALL FFTC(B1(1, 1) ,B1(2,1) ,LN,2) 

DO J ■ 1 , FRAME_SIZE 

ROLLER (1, J) -r ROLLER ('1, J) + AREAL(C1(J)) 
ROLLER(2,J) * ROLLER(2,J) + AIMAG(C1(J)) 
END DO 
END DO 
END DO 

MOD MAX - 0 

DO I ■ 1 7 FRAME_SIZE 

ROLLER(l,I) - ROLLER(l,IJ / (NEX_IN_RAW * 2) 
ROLLER(2,I) - ROLLER(2,I) / ( NEX_IN_RAW * 2) 

ROLLER (3,1) - SQRT(ROLLER(1 > I)*ROLLER(1,I)+ROLLER(2 / I)*ROLLER(2 



IF (ROLLER(3,I) .GT. MOD_MAX) THEN 

MOD MAX - ROLLER(3,I) 
END IF 

IF ( (ROLLER(l p I) .EQ. 0) .AND. (ROLLER(2,I) * EQ. 0)) THEN 

ROLLER (4 , 1 ) - 0 
E LS E 

ROLLER(4,I) • ATAN2 (ROLLER (2,1) , ROLLER (1,1)) 
END IF 
END DO 

MOD_25 - MOD_MAX * .25 

start roll - 1 
stop Roll • -42 

do while ( (start roll . lt. frame_size - 16) « and. 

# (stop_roll .lt. 0)) 

do while ((roller (3, start roll) .lt, mod_25) .and . 

# (start roll • lt* frame_size-16) ) 

START_ROLL - STA5T_R0LL + 1 
END DO 

IF (START ROLL .LT. FRAME_SIZE - 16) THEN 

IF ( (ROLLER (3, START ROLL+5) • GE . M0D_25) .AND. 

# ( ROLLER ( 3, 5TART~ROLL+6) .GE. MOD_25) .AND. 

# (ROLLER ( 3 , START~ROLL+7 ) .GE. MOD_25) • AND. 

# (ROLLER (3 , START_ROLL+8 ) . GE. MOD_25) .AND. 

# ( ROLLER ( 3 , START ROLL+9) .GE. MOD_25) .AND. 

# ( ROLLER ( 3 , START~RO LL-f 1 0 ) . GE. MOD_25) .AND. 
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I / C «S^^f ^ ' f — * > « GE * «°D 25) .AND, 

t ^^Ste'S****-* 01 ** 1 *) .GE. MOD-25) .AND 

t (ROLL£R( 3 , START_ROLL*hl 3 ) . GE. HOD~25> !and 

* ( ROLLER ( 3 f START ROLL+ 14 ) .GE. MOD^25> I AND 
I (ROIiER(3,STARTlR0LL + a5 . GE. MOD-fl SnSI 

* (ROLLER (4, START ROLL+5) ,LT. PI9)~ .AND 
I <ROL^R(4,START:ROLL-k6 .LT. Pli j&S: 

! / ( IS^S / (4 ' S7ART - R0I ' L+7 > * LT - P « -AND. 

* (ROLLER (4, START ~ROLL+8) .LT. PI9) .AND 

* ROLr£R( 4> STARTlROLL + 9> .LT. PI9 . ^51 

* JS8H5^fE& MLI * 10) - LT - P "> -AND. 

* <KKSH5 rSMOT - R0LL * I1 > «LT. PI9) .AND, 
t CROIXER(4rSTART.ROLL+12> . LT. PX9> .AND. 

* ( ROLLER (4 > START_ROLL+ 13 ) .LT. PI9> .ANIV 

* (ROLLER ( 4 £ START ROLL* 14 ) .LT. PI9 AND 

* OT ^S^?, (4 ' START I R0LL+15 PI ) THEN 
START_ROLL - STAR? ROLL + 5 - 

„^SI OP - ROLL * START ROLL + 10 
ELSE ~ 

START_ROLL - START ROLL + 15 
END. IF ~ 
END IF 
END DO 

IF (STOPJROLL .GT. 0) THEN 
Exxyi *.0.0 
Exi m 0.0: 
Eyi » o.O 
Exl2 m (KG 

DO J ■ START ROLL r STOP /ROLL 

Eyi « Eyi +- ROLL£R(4,J> 
Exi2 « Exi2 + J * J* 
END DO 

1F K&Tt Z Zi2^ (Exi * ? xi > - NE - °> THEN 

EISB - " (li*Exiyi - Exi*Eyi) / (ll*Exi2 - Exi*Exi) 

ANGLE ROLL w "999999999 
. END IF ~ 
END IF 
END IF i DO_ROLL 

IF (DOJSHIFT* THEN 

IF (IE^!NE?0) ( ^ TEP * P ^ E - FRAME -. FR °M^TEP i BUFFER, ICHAN1 , IERR) 
^PRlll^*"^ " X> * AND * (IERR .NE. 24) ) THEN 
*PHASE_FRAME FR^ c6de! ^IrI^"^ ^ BLOCK NUMBER: ',23-3TE? 

Else 

print*, • 1 . 

PRINT*, 'ERROR, END OF FILE. 1 
END IF 
RETURN 
END IF 

B_C « FOFF * ((NAV_ECHO-l) + PKASE_NEX_FROM * NECHOES) 



WO 90/05920 



PCT/US89/04945 



49 



1999 



2000 



IF (MOD2(PHASE_NEX_FROM,TWO) .EQ. 1) THEN 

CHOPPER FAC - -1 
ELSE 

CHOPPER FAC - X 
END IF 

L * FRAME SIZE/2 + 1 

DO J - 1+1 C , FRAME SIZE+B C , 2 

C1(L) - 5MPLX( BUFFER (J )"* CHOPPER FAC r 
BUFFER (J+l) * CHOPPER~FAC) 

L - L + 1 
END DO 

L - 1 

DO J » FRAME SIZE+l+B C , FRAME SIZE*2+B C , 2 
CI(L) - CMPLX ( BUFFER ( J ) * CROPPER FA£, 
BUFFER (J+l) * CHOPPER~FAC) 

L « L + 1 
END DO 

CALL FFTC(B1<1,1) ,B1(2.1) ,LN,2) 

MOO MAX - 0 

DO J - 1 , FRAME SIZE 

ROLLER(l , J) • AR£AL(C1(J)) 

ROLLER(2,J) • AIMAG(CKJ)) 

ROLLER(3,J) - SQRT (ROLLER ( X , J ) * ROLLER ( I # J ) + ROLLER ( 2 , J ) * ROLLER ( 2 

IF (ROLL£R(3„J) .GT. MOD MAX) THEN 

MOD MAX - ROLLER(3,J) " 
END IP 

IF ((ROLLER(X,J> .EQ. 0) .AND- (ROLLER (2, J) . EQ. 0)) THEN 

ROLLER(4,J) • 0 
ELSE 

ROLLER(4,J) - AT AN 2 (ROLLER (2 , J) , ROLLER ( 1 , J) ) 
END IF 

END DO 

MOD_10 - MODJKAX * . 10 
WRITE (6, 1999) 

FORMAT (/, 1 HERE ARE PHASE ANGLES CLOSE TO 0 1 ,/) 

DO J - 1 , FRAME SIZE 

IF ((ROLLER (3, 7) * GE. MOD_10) .AND. 
(ROLLER(4,J) .GT. -0.1) .AND. 
(ROLLER (4 , J) .LT. 0.1)) THEN 
WRITE (6, 2000) J, ROLLER ( 4 , J) 

FORMAT ( 1 COL - M3,' ANGLE - »,F6.4> 

END IF 
END DO 



WRITE(6,2020) 
PHASE COL - -4 2 
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"J^MS^** * LT - « -° R; FRAME_SIZE>> 

3000 FOKMATC'+ENTER THE COLUMN TO USE FOR THE PHASE. SHIFTS <GT» , $ 

^j^SS-COL - INT2 (GET_IO ) 

PHAS£_R£F - ROLLER (4, FHASE_COL) 
END IF 1 DO_SHIFT 
IF (DO_DELTAS) THEN 

XFOT?S?S) ( ^ TEP#FRWa Lrao«' STEP , BUFFER r ICHANI ; I ERR) 
IF pRINT^r^ E ' -AND, (I ERR .NE. 24)) THEN 

E : ' , IERR PRINT ,ERHOR DEFECTED WHEN READING BLOCK NUMBER: M, - CC3 

ELSE * 

PRINT*,' • 

^J?***!**' 'ERROR, END OF FILE. • 
E ND I F 
RETURN • 
END IF 

B^C - FOFF * ( CNAVjECHO-1) + NEX_FROK * NECHOES) 

XT «SppMl^™ 0 ?£ TWO) - EQi l > THEN 
ELSE 

CHOPPER FAC * 1 
END IF - 

. L/ - FRAME SIZE/2 + 1 . 

D° * • 1+B C , FRAME SIZE+B C , 2 

CI(L) - &!PLX(BUFFER(,r )~* CHOPPER FAC, 

* » _ BUFFERCJ+I> * CHOPPER"FAC> 
** * t + I 

END DO 

Ir * 1 

D °cir^ F 2 A S5w^?S?^ff r ^' FRAME SIZE*2+BC , 2 
CI CD • CMPLXCBUFFERCT ) * CHOPPER FAC r 

* . BUFFER(Ji-l) * CHOPPER~FAC) 
U ■ m u + X "~ 

END DO 

CALL FFTC(BX(X,1) ,BIC2 r L>,IN,2) 
Ir - 1 

DO J » 1 / FRAME SIZE/2 
RR * AREAL(C1(J) ) 
RI - AIMAG(CICJ)) 

VIEWXCL) - SQRT (RR*RR+RX*Ri) * ZBUF (L+HIGH HALF) 
It » L. * 1 — 
END DO 

L * LOW_HALF 
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DO J - FRAME SIZE/2 + X , FRAME_SIZE 
RR - AREALlCl(J)) 
RI - AIMAG(C1(J*)) 

VIEWl(L) - SQRT (RR*RR+RI *RI ) * ZBUF (L+HIGH_HALF) 
L ■ L + 1 
END DO 
END IF ! DO_DELTAS 

WRITE(6,2020) 

IF ((DO DELTAS) .OR. (DO SHIFT)) THEN 
I - 2l 
K - 0 
IERR - 0 

DO WHILE (I • LT. STOP_BLOCK) 
K • K * 1 

CALL LIBSRDB ( I , STEP, BUFFER, I CHAN1 , IERR) 
IF(IERR.NE.O) THEN 

IF ((IERR .NE. -1) .AND. (IERR . NE. 24)) THEN 
PRINT*,' 1 

PRINT *, 'ERROR DETECTED WHEN READING BLOCK NUMBER: '.I. 

ODE: 1 , IERR 

ELSE 

PRINT* , 1 • 

PRINT*, 'ERROR, END OF FILE. 1 
END IF 
RETURN 
END IF 

WRITE (6, 2010) K 
2010 FORMAT ( ' <GT> 1 ,13,$) 

IF (DO DELTAS) THEN 

DO LT - 0 , NEX IN RAW - 1 

B_C « FOFF * ( (NAV_ECHO-l) + L2 * NECHOES) 

IF (MOD4(L2,2) .EQ. 1) THEN 

CHOPPER FAC - -1 
ELSE 

CHOPPER FAC - 1 
END IF ~ 

L - FRAME SIZE / 2 * 1 

DO J - 1*5 C , FRAM£_SIZE*BJZ , 2 

C1(L) - CMPLX( BUFFER (J ) * CHOPPER_FAC , 
f BUFFERED * CHOPPER_FAC) 

L - L + 1 
END DO 

L « 1 

DO J - FRAME SIZE+l-B C , FRAME SIZE* 2 -B_C , 2 
C1(L) « CMPLX ( BUFFER (J- ) * cEOFPER_FAC , 
4 BUFFER (J+l) * CHOFF£R_FAC) 

L » L + 1 
END DO 

CALL FFTC(B1(1,1) ,B1(2,1) ,LN,2) 
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IF (DO SHIFT) .THEN 
PKASfC - PHASEC «• 1 

IF ( (AR£AL(Cl( PHASE COL) ) . EQ. 0) .AND. 
(AIMAG(Cl(FHASEj:OL)) .EQ. 0)) THEK 
PHASES (PHASEC) • 0 - PHASE REF 
ELSE 

PHASES (PHASEC) - 

ATAN2 (AIMAC (CI (PHASE. COL) ) ,AREAL( CI (PHASE COL))) - 
# PHASE REF " ~ 

END IF ~ 
END- IF I DO_SHIFT 

L * 1 

DO J • i , FRAME SIZE / 2 
RR » AREAL(CI C3F)") 
RI ■ AIHAG CCl(JT) ) 
VIEWN(L) ~.SQRT(RR*RR+RI*RI) 

END DO 

L » LOW HALF 

DO J - ERAKEJSIZE / 2*1, FRAME SIZE 

RR - AREAL(C1(J>) 

RI AIMAG(C1(J>> 

VIEWN(L) » SQRT (RR*RR+RI*RI) 

L « L «■ I 
END DO 

VECTORC - VECTORC * 1 

DO L - LOW CON , HIGH CON 

DO M - L5«_HALF , HIGH HALF 

THE KODS(M) ~ VXEttl(R+t,> * VIEWN(M) 

MOD _MEANS(L) - 0.0 

• DO K « LOW HALF , HIGH HALF 

MOOJMEANSXL) - MOD MESIAL) * THE HODS(M) 
£NO DO " ~ " 

MOD_MEANS(L) « MO DEMEANS (L) / 256 ! (HIGH.CON - LOWJTCN - 
END DO 



MAXOV1 - -42 

DO L" LOW CON , HIGH CON 

if (mod Means (L) .gt; maxovi) then 

MAXOVT * MOD MEANS (L) 
DELTAS (VECTORC) = L 
END IF 
END DO 
END DO i FOR EACH NEK 
END IF I DO DELTAS 
I - I + STEP 
END DO 

END IF ! DO DELTAS OR DO SHIFT 
2020 FORMAT ( * •) ~ 

WRITE (6, 2020) 

RETURN 
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END 



END OF DO MARCO. F77 SOURCE FILE 



Those are the two main components of the MARCO program. The following 
two are the main components of the POLO program, the program that 
actually changes the raw data to implement the motion crrections 
mentioned above. The main POLO program gets some information from 
the user and reads the text file created in the MARCO program, and the 
routine DO_POLO implements the changes. 



BEGINNING OF POLO.F77 SOURCE FILE 



C 
C 



COPYRIGHT - MAYO FOUNDATION FOR MEDICAL EDUCATION AND RESEARCH - 1983 
POLO 

THIS PROGRAM TAKES THE DATA GENERATED BY MARCO AND APPLIES 
THEM TO A RAW FILE, WRITTING THE NEW DATA TO A RAW FILE 



IMPLICIT NONE 
CHARACTER* 12 8 SOURCE FILE 
CHARACTER* 12 8 ASCII_2ATH 
CHARACTER* 12 8 D£SP_PATH 

INTEGER* 4 IERR 
INTEGER*4 ICHAN1 
INTEGER* 4 ICHAN2 



Name of source file 

Name of file of ascii values 

Name of file with displacement 

correction 

Error flag 

Channel number for source file 
Channel for displacement correction 



c. 
C. 
C. 
c. 



c. 
c. 



The following are identifiers that are used with the routine 
GET — COLLECTION , which, gets some information about the raw data. 

CHOPPER ! If the file has chopper data 

HOAP ! If the data has been processed by 

! array processer 
CHOPLETSON .! If the file has chopletson data 

SHIM, HNEX , STRIP, D3FFT, GOT_IT 
EXORCIST , HI_LO^SORT, NO_SORT 
SKIP920 

BASELINE SKIP i Number of blocks of baseline data 
FRAMES PER SLICE I Frames per slice . 
DATA BlOCKS* I Number of data blocks 
HAVS~ i Number of averages 

NECHOES I Number of echoes 

FRAME SIZE 1 Frame size 



LOGICAL* 2 
LOGICAL*2 

LOGICAL* 2 
LOGICAL*2 
LOGICAL*2 
INTEGER* 2 
INTEGER* 2 
INTEGER* 2 
INTEGER*2 
INTEGER* 2 
INTEGER* 2 
INTEGER* 2 



Number of bytes in opened file 



INTEGER*4 I BYTES 

These next two relate to the ascii file that holds information 
from the marco program. If the file was written in an older 
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c. 
c. 

c. 



c 
c 
c 



c. 
c 



c 
c 

c 
c 

c 

c 

c 

c 



142 



version of the program, only delta values and frame numbers are 
line h 2e5i£e w£h "iKFoT"^ 0th * r inf0 ™*"°n, **d the first 

LOGICAL*2 HEW VERSION 
CHARACTER*5 lS THIS INFO 



INTEGER* 2. VECTORC 
INTEGER* 2 NEX 
INTEGER* 2 NEX IN RAW 



The near two help make the 
INTEGER*2 STARTO 
CHARACTER UNIQUE 

INTEGER** IXX r I, K r J 



! Counter into VECTOR, below 

i Number of excitations 

> Number of excite ions in the raw 

! file. Data that have been processed 

J by the array processer has only one 

1 while other chopper data have 2 or 4 

new raw data file's name. 



! counters 



Th 



e next fifteen store information from the ascii data file. 



INTEGER* 4 DELTAS CX024) 

REAL PHASES (1024) 
LOGlCAL*2 DO DELTAS 
LOGICAL*! DO - ROLL 
LOGICAL* 2 DO~SHXFT 
INTEGER* 2 FRjME_DELTA 

INTEGER* 2 START_DELTA 
INTEGER* 2 STOP^ DELTA 



INTEGER* 2 START_ROLL 

INTEGER*2 STOP_ROLL 

INTEGER* 2 FHASE_COL 

INTEGER* 2 P HAS ENFRAME 

INTEGER*2 IMAGE ECHO 
INTEGER* 2 NAV ECHO 
REAL ANGLE_ROEL. - 

REAL* 8 DELTA FACTOR 
REAL* 8 PHASE~FACTOR 
INTEGER* 4 PAESE^F 

WRITE(6,142) 
FORMAT C ' 1 * # / /////) 

rxi - o 

GOT IT * .FALSE. 

do While ((.not. got m 



The delta values r how far to shift the 
data 

The phase shift in a selected column 

To shift by the delta values 

To do the roll correction 

To do phase shift correction 

What frame was used as a reference for 

getting the deltas values 

The sample number at which the 

determination of the delta values 

began 

The sample number at which the 
determination of the delta values 
ended 

Column at; which roll angle 
determination began 
Column at which roll angle 
determination ended 
Column used for phase shift 
determination 

Frame used for phase shift 

determination 

which echo is the image 

Which echo is the navigator 

The roll angle 



•AND. (Ill .LT. 10}) 
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GOT IT - .TRUE. 
Ill - - III + 1 
WRITE(6,3010) 

3010 FORMAT (•-•-ENTER SOURCE FILE 

«5T>' ,$) 

R£AD(*, • (A) • ) SOURCE_FILE 
DO I- 1 , 128 

IF(SOURCE FILE(i:I) .EQ.«<40>") SOURCE FILE (I : I) «"<0>" 
END DO 

CALL LIBSGOPEN (SOURCE FILE , 0 , ICKAN1 , IBYTES , IERR) 
IF (IERR.NE. 0) THEN 
IF (IERR. EQ. 21) THEN 

PRINT *, "SOURCE FILE DOES NOT EXIST." 
PRINT* , CHAR ( 7 ) 
ELSE IF (IERR . EQ. 20) THEN 

PRINT*, 1 ILLEGAL FILENAME CHARACTER. • 
PRINT* , CHAR ( 7 ) 
ELSE IF (IERR .EQ. 172) THEN ! ZERO LENGTH FILENAME 

CONTINUE 
ELSE 

PRINT*, 'ERROR IN OPENING SOURCE FILE, CODE: 1 , IERR 
PRINT* , CHAR ( 7 ) 
END IF 

GOT IT • .FALSE. 
END if 
END DO 

IF (III .GE. 10) GO TO 9000 

C. . . 

C... THIS GETS SOME INFORMATION ABOUT THE RAW DATA FILE BY READING 
C... PARTS OF THE FILE HEADER. 

CALL GET COLLECTION ( I CHAN 1 , CHOPPER, SKIM/ HNEX , STRIP , D3 FFT , 

# EXORCIST , HI LO SORT, CHOP LETS ON, NOAP, NO SORT, 

# BASELINE SKTP , 3KIP9 2 0 , FRAMES PER_SLICE, DATA_BLOCKS , 

# FRAME SlZE,NECHOES,NAVS,IERRl 
IF (IERR . NE."~0) T H EN 

PRINT*, 'ERROR FROM GETTING COLLECTION TYPE, CODE : * , IERR 
GO TO 9000 
END IF 

IF ((CHOPPER) .AND. (.NOT. NOAP) ) THEN 

NEX IN RAW - 1 
ELSE TF 7 (CHOPPER) .AND. (NOAP)) THEN 

NEX IN RAW - 2 * NAVS 
END if 

III - 0 

GOT IT - .FALSE. 

do While ((.not. got it) .and. (in ,lt. ion 
in - in + i 

WRITE(6,3015) 

3 015 FORMAT ( •♦ENTER FILE OF DELTA VALUES* 

<GT>« ,$) 

READ(*, • (A) ».) ASCII_PATH 

DO I- 1 , 128 
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^FXASCII.PATHCf:!, . EQ . , <40>11> ASCII_PATHCI:t } -< 0 >» 

^^ I !^?i ) II ^ Clr - P * T «^*WS-.OU,.,IO S TAT.XE R R, 
e J|INT *, 'ERROR DETECTED DURING PILE OPEN. CODE : • ZtBR 

END DO 

IF (III .GE. 10) GO TO #000 

aafcfssfia^ SLICE . NAVS . 2 

. C... HAS MORE INFORMATION . NO INFORMATION. THE NEW VERSION 

4241 F^a^aI> 41>ZS - TH1S - IKFO 

enS^F™ 101 ' " ' FALSE - 

REWIND(42) • 

IF (NEWVERSION) THEN 

END DO 

EL ^ iLJ5 XS IS OLD VERSION 

DO_ROLL - . FALSE \ 
?O r fHIFT - « FALSE. 
S TXRT DELTA » -1 

stopJJelta • -1 

00 •£ " 1 > VECTORC 

END DO 
END IF 

CLOSE (42, IOSTAT-IERR) 
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IF ( D0_ DELTAS ) THEN 
UNIQUE - •** 

DO WHILE ( (UNIQUE . NE. 'Y') .AND. (UNIQUE . NE. 1 N 1 ) , AND. 

* „ (UNIQUE .NE. «y») .AND. (UNIQUE .NE. 'n')) 

WRITE (€, 4000) 

4000 FORMAT ( *+DO YOU WANT TO DO THE DELTA SHIFT CORRECTION ? (Y/N1-- 

— <GT>',$) 

READ ( 5 , 3999 ) UNIQUE 
3999 FORMAT (A) 

END DO 

IF ((UNIQUE .EQ. *N • ) .OR. (UNIQUE .EQ. 'n')) THEN 

DO DELTAS - .FALSE. 
END TF 

IF (DO DELTAS) THEN 
WRITE(6 f 4500) 

4500 FORMAT ( ■ + ENTER THE DELTA SHIFT FACTOR 

<GT> ' , $ ) 

K - PARS£_F ( 1 , DELTA FACTOR) 
IF (K .EQ. 0) DELTA~FACTOR - 1 
END IF 
END IF 

IF (DO ROLL) THEN 
UNIQOE - 

DO WHILE ((UNIQUE .NE. • Y ' ) .AND. (UNIQUE .NE. «N») .AND. 

# (UNIQUE .NE. 'y*) .AND. (UNIQUE .NE. , n")) 
WRITE(6,4001) 

4001 FORMAT ( '+DO YOU WANT TO DO THE ROLL BACK CORRECTION ? (Y/N) 

<GT>',$) 

READ (5,3999) UNI QUE 
END DO 

IF ((UNIQUE .EQ. »N*) .OR. (UNIQUE .EQ. I n , >) THEN 

DO ROLL - .FALSE. 
END IF 
END IF 

IF (DO SHIFT) THEN 

uniqOe - 

do while x (unique .ne. 'y') .and. (unique .ne. 'n») .and. 

# (UNIQUE .NE. •y') .AND. (UNIQUE .NE. 'n')) 
WRITE ( 6 , 4 002 ) 

4002 FORMAT < *+DO YOU WANT TO DO PHASE SHIFT CORRECTION ? (Y/Ni 

~<GT> • , 5 ) 

READ (5, 3 999) UNIQUE 
END DO 

IF ( (UNIQUE .EQ. 'N') -OR. (UNIQUE .EQ. 'n')> THEN 

DO_SHIFT - .FALSE. 
END IF 

IF (DO SHIFT) THEN 
WRITE(6,4510) 

4 510 FORMAT ( ' + ENTER THE PHASE SHIFT FACTOR 

<GT>' ,5) 

K - PARSE F(l , PHASE FACTOR) 
IF (K .EQ? 0) PHAS E~FACTOR - 1 
END IF 
END IF 
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.IF (CHOPLETSOK) THEN . 

HEX ■ 1 
ELSE 

KEX - NAVS * 2 
END IF 

IF (HNEX) THEN 

J^55;*^ TaiAB£E T0 PROCESS HALF NEX DATA. 1 
GO TO 9000 
END XF 

IF «X^Ji2T* . ^ (CHOPpER > .OR. (CHOPLETSON) . 0 R. (NOAP) V ) THEN 

GO TO 9 000^ CH0PPER «. N0AP < OR CHOPLETSON, UNABLE TO PROCESS. 1 
END IF 

I - 1 

DO^WHILiEl^CCI -LT. X22) .AND. (S0URCE_FILE(I : I) .NE. CHAR(O) ) J 
END DO 

IF (I .LT. 122) THEN 

IF (CSOURCE_FILE(T-4:I-l> <EQ. '.RAW') .OR. 
(SOURCE FILE (1-4: 1-1* . EQ . ».raw-)) THEN 
IF (<SOUR^E.FILE(I-5:I-5) . LT. *0») .OR. 

e J ISH* 0 ^" 1 * 51 (1-5:1-5) .GT. '9*)) THEN 
STARTO m X-5 
ELSE 

STARTO m 1-4 
END IF 
ELSE 

STARTO -X 
END IF 

UNIQUE 'm. 

DO WHILE . ( (UNIQUE . LT. «A») .OR. 

I C JHS5SH5 -SI* ' 2,) • AND - (UNIQUE . LT. 'a')) .OR. 

* (UNIQUE .GT. ■Z 1 )) 

WRXTE(6,3019) . 
3019 F0R1«lT(«+ENTER CHARCTER TO UNIQUELY IDENTIFY THE OUTPUT FX Li 

—————— <GT>' r S) 

RSADC*, • (A) l y UNIQUE 
END DO 

^P^PATH^-^OURCE_FILE»( 1 : STARTO-1 ) //UNIQUE// 1 . DSP 1 

DESF_PATH(I:I) '« CHAR(O) 
END DO 
ELSE 

PRINT*, 'ERROR, SOURCE FILENAME IS TOO LONG. • 
GO TO 9000 
END IF ■ 

C. . . THIS MAKES A NEW FILE FOR THE RAW DATA. 

CALL MAKE NEW UDF ( DESP PATH, IERR) 
IF (IERR .NE. 0) THEN ~" 

SF2 T ** ,£RROR DETECTED DURING DESP FILE MAKE/ CODE : ' r IERR 
GO TO 9000 
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ENO IF 

CALL LI B$ CO PEN (OESP PATH, 0 , ICHAN2 , IBYTES , I ERR) 
IF ( I ERR « NE . 0) THE? 

PRINT *, "ERROR DETECTED DURING DESP FILE OPEN, CODE : MERR 

GO TO 9000 
END IF 

IF (CHOPPER .AND. NOAP) THEN 

WRITE (6, 3040) 
ELSE IF (CHOPPER .AND. .NOT. NOAP) THEN 

WRITE (6, 3050} 
END IF 

WRITE (6 , 3060) FRAMES PER SLICE, FRAME SIZE 
WRITE ( 6 , 3 07 0 ) NECHOES , NEX 

3 040 FORMAT (• SOURCE FILE: CHOPPER AND NOAP ' ) 

3050 FORMAT (• SOURCE FILE: CHOPPER AND AP«) 

3060 FORMAT ( 9 FRAMES PER SLICE - • r I3 ,T32 , 1 FRAME SIZE - M3) 

3070 FORMAT ( 9 ECHOES - * , I3,T32, 'NEX - ' ,13,/) 

C... DO_POLO IS THE ROUTINE WHERE ALL THE NUMBER WORK IS DONE 

CALL DO POLO (DELTAS, PHASES, DO DELTAS, DO ROLL, DO SHIFT, 
« START DELTA, STOP DELTA, ANGLE ROLL, NEX IN RAW , NECHOES , 

# FRAME'S XZE , FRAME? PER SLICE, BAT A BLOCKS, IMAGE ECHO, 

# ICHANl r I CHAN 2 , DELTA FACTOR , PHASE~F ACTOR , IERR) ~ 
IF (IERR . NE. 0) THEN 

PRINT*, 'ERROR FROM DO POLO, CODE-", IERR 
END IF 

9000 CONTINUE ! GONE TO FOR ERRORS 
END 



END OF POL0.F77 SOURCE FILE 



BEGINNING OF DO POL0.F77 SOURCE FILE 



COPYRIGHT - MA^O. FOUNDATION FOR MEDICAL EDUCATION AND RESEARCH 

SUBROUTINE DO POLO (DELTAS , PHASES , DO DELTAS , DO ROLL , DO_5HIFT , 

# START - DELTA , STOP DELTA, ANGLE ROLL, NEX IN RAW, NECHOES , 

# FRAME'S IZE, FRAMES PER SLICE, DATA_BLOCKS ,TMAGE_ECHO , 

# ICHANl, ICHAN2. D£LTA_FXCTOR, PHASE FACTOR, IERR) 
IMPLICIT NONE 

INTEGER*4 IERR, DELTAS (*), B_C 
REAL PHASES (*) 
INTEGER* 2 IMAGE ECHO 
INTEGER* 2 BUFFER ( 32768 ) , IER 

INTEGER* 4 STEP, I, J, K r L, L2 , ICHANl, ICHAN2 
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c. 
c. 
c, 
c. 



c. 
c. 



SSSR SS&^S^*' "MBW, STOP.BXOCK, FRAME SIZE 
!ara«L S I$g- DELTA ' ST0P ^ DE "*' NEX.INRAW 



REAL*8 DELTA FACTOR 
REAL** PHASE^FACTOR 

INTEGER* 4 LN, GET_I - 

L0GICAL*2 DO_DELTAS, DO.ROLL, DO_SHIFT 
REAL RR, RX, NEW, SNEWV CNEW 

COMPLEX CI (512) 
REAL BI(2,512) - 
EQUIVALENCE [Cl # B%) 

COMPLEX C2(512) 
REAL B2 C2 r 512) 
EQUIVALENCE (C2, B2 j 

INTEGERS STARTJkT, END^AT 
WRITE C 6, 2 020) 

START AT m i 
END_AT • FRAME_SI2E 



|: : : SShSa^SffiS. C0PIED *— 152 ? OORCS FItE » ™= 

raiOT*. 'ERROR FROM MOVING HEADER, CODE r IER 
END IF 

X2SS-i S *ffP_OFTSET THAT IS USED WHEN DATA ARE RFAD rami » 
B0TFER. THERE ARE 2 a-HRfe INTEGERS FOR EACH SAMPLE™ 

FOFF » FRAKE_SIZE * 2 

f??e I S.2? E ^5! BER - OP 512 BWE BLOCKS TO READ FOR EACH VIEW 
LN IS USED WITH THE FFT ROUTINE, 2**LN ■ FRAKE/.SIZE 

^inS** 8 - 8 ? 28 - E Q- 5 ") THEN 
STEP - 4 * NECHOES * HEX_IN_RAW 

Bl fSr£ F ffRAKE-SIZE -EQ. 256 > THEN 
STEP - 2 * NECHOES * HEX IN RAW 

Ur ^8 " 

ELSE " 

IERR « 4201 
RETURN 
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END IF 

C. . . 

C... STOP BLOCK IS THE LAST BLOCK THAT SHOULD BE READ. 
C. . . 

STOP_BLOCK - 28 * STEP * FRAMES_PER_SLICE 

IF (STEP * FRAMES_PER SLICE .GT. DATA BLOCKS) THEN 

PRINT*, 'ERROR, NOT ENOUGH DATA BLOCKS* 1 

I ERR - 4200 

RETURN 
END IF 

WRITE(6,2020) 

C. . . 

C... I IS THE BLOCK TO BE READ, IT IS INCREMENTED BY STEP BELOW. 

C. . . K IS THE VIEW NUMBER CURRENTLY BEING DEALT WITH. 
C. . . 

1-28 

K - 0 

I ERR - 0 

C. . . THIS LOOP IS DONE ONCE FOR EACH VIEW 

DO WHILE (I . LT. STOP BLOCK) 
K - K + 1 ~ 
WRITE (6, 2010) K 
2010 FORMAT ( 1 *+— <GT> ' , 13 , $ ) 

CALL LIBSRDB ( I , STEP , BUFFER , I CHAN 1 , IERR) 
IF(IERR.NE.O) THEN 

IF ( (IERR .NE. -1> .AND. (IERR .HE. 24)) THEN 
PRINT*," 1 

PRINT *, 'ERROR DETECTED WHEN READING BLOCK NUMBER: M. • CCD 

E:',IERR 

ELSE 

PRINT*, ■ • 

PRINT*, 'ERROR, END OF FILE. 1 
END IF 
RETURN 
END IF 

C. . . 

C. THIS LOOP IS DONE FOR EACH NEX IN THE RAW DATA. CHOPPER- AP 

C. . . DATA HAVE ONE NEX IN THE RAW DATA, CHOPPER NOAP DATA HAVE 

C... TWO OR FOUR, OR MAYBE MORE. 

C. . . 

DO L2 - 0 , NEX_IN_RAW - X 

C. THE DELTAS ARRAY HAS THE AMOUNT OF SHIFT THAT IS NEEDED. 

IF ( DELTAS ( K*NEX_IN_RAW- 1+L2 ) .NE. 0) THEN 

C... B C IS A BUFFER COUNTER THAT IS HOW FAR INTO THE BUFFER 

C. . . T5 BEGIN TO READ FOR THE NEX 

B_C - ((IMAGE_ECHO-l) + L2 * NECHOES) * FOFF 
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c. . 
c. . . 



SAMPLE 129 GOES TO wMtiS?*! T ?,S2 S I TI0,r AND 
g- :: ARRA* HAS 2 Z^raS^NTECER^* **** «*"».» THE BUFFER 

E --FRAME SIZE/2 "♦ I 

,»f,r. 1 "'" 5 - C » FRAME SIZE+B C 2 

Cia). ^^UCfWFFERwf.BBFFER?^!,, 
END DO 

END DO 

CALL FFTCCBl(l,lX r Bl(2 f l),LM,2J 
xr r C -NOT* DO_DELTAS) THEN .. 
ELSE 

ETO J p f 0ELTA S^*NEX.IK.RAW^ tt ] * DELTA_TACTOR) 1 X 
J » STAHT^AT * 1 

DO^WHILE^aj ,LT. END.AT) .AND. (L .LT. 0) ) 

en§*&*^- siz& * } " «W 

l V5 1 £Vf r ' LT - ^-^J -AKD' (L .LT. FRAME SIZE) ) 
J - J + 1 

DO WHILE (J .LT, END AT) 

^CZa-FRAME.SIZE) . CI (J) 

IF (CDO_ROLE) .OR. fDO_SHlFT)) THEN 
'. NEW *'*Q 

IF (DO_ROLL) THEN 
% NEW » NEW. ♦ (BSLTAS(K.NE^IN_RAW-l + L2) • DELTA FACTOR ) 

END IF ANGLE^ROLL 

IF (DO_SHIFT) THEN 

NEW a HEW + (PHASES ( K*NEX_IN.RAW-1 + L2) * EHASE.FACTCR) 
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: ' , IERR 



END IF 

SNEW • SIN (HEW) 
CHEW m COS (NEW) 

DO J m X , FRAME SIZE 
RR - B2(1,J) " 
RI - B2(2,J) . 

B2(1,J) « RR * CNEW - RI * SKEW 
B2(2,J) - RR * SHEW + RI * CNEW 
END DO 

END IF J DO_ROLL OR DO_SHIFT 

CALL FFTC(B2(1,1) ,B2(2,1) ,LN,-2) 

L - FRAME SIZE+l+B C 
DO J - X 7 FRAME SIZE/2 

BUFFER (L ) - lNT2(AR£AL(C2(J) )) 
BUFFER (L+l) - INT2(AIMAG(C2(J))> 
L - L + 2 
END DO 

L - 1+B C 

do j - Frame size/2+i , frame size 

BUFFER (L T - INT2(AREAL(C2?J))) 
BUFFER(L+1) - INT2(AIMAG(C2(J))) 
L - L + 2 
END DO 
END IF 

END DO ! FOR EACH NEX IN RAW 
CALL LIBSWRB (X,STEP,BUF?ER, ICHAN2 , IERR) 
IF (IERR «NE. 0) THEN 
PRINT* , ' ' 

PRINT *, 'ERROR DETECTED WHEN WRITTING BLOCK NUMBER: »,!,• CCDS 

RETURN 
END IF 

r - I + STEP 
END DO 



2020 FORMAT ( * •) 

WRITE (6 ,2020) 

RETURN 
END 



END OF DO P0L0.F77 SOURCE FILE 



That concludes th« four major components of the MARCO and POLO pro- 
grams* Below are other source file used in one, some, or all of the 
above source file. 



BEGINNING OF PARCE.F77 SOURCE CODE 
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j COPYRIGHT - KAYO FOUNDATION FOR MEDICAL EDUCATION AND RESEARCH 198* 
FUNCTION PARSE I (HOW, NUHS) 

PARS! I r HOW r NUMSf*> 
CHARACTERS 0 ST£ 1 9 

IOGICAL*! NEGATIVE r DONE 

READ (5, 1000) STR 
1000 FORKAT |A8 0 ) 

K - 0 
I » 1 

DONE • .FALSE* - 

DO j mtt - (( f H0T, 0QR8) .AND* (K*LT. HOW)) 

DO WHILE ffl .LT. 80) . AND, m(Zz I) .NE* .XNQ 

* crsnux-x! lt -on (s £ (rs ,£l»;?^ ,+t > -and? 

I - X-*™ * -OR. .(ST»(X*XJ .GT* *9*))) 

EKD DO 

. IF (I .LT* EO) THEN 

IF .J STR < l5 *> *EQ* THEN 
NEGATIVE .TRUE. 
I » I + 1 

BI SLJjjK TR ' Jtr:I j- -EQ. THEN 
. . NEGATIVE - . FALSE*. 

1*1 ♦ 1 

ELSE 

^NEGATIVE - .FAISE. 

^k"*?^^ ,0,) M * fSTRCIrl) . LE. . 9 ')) THEN 

DO WHILE fJ^; x ~> rf » 5i J«CX.Xl MM, .AND* 

X Z J * 1° * ^ ICHAR ( S TR(X-X)) - 48) 
END DO 

ELSE 

DONE - .TRUE. 

END zr 

IF (NEGATIVE)' THEN 
^HUMSOC) - J . f-i, 

ELSE 

DONE • '.TRUE. 

END IF 
END DO 
PARSE I - K 
RETURR 
END 
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FUNCTION PARSE F( HOW, HUMS) 
INTEGER* 4 PARSE F, HOW. K. J 
REAL*! HUMS (*) ~ ' 
REAL* 8 J2 
CHARACTER* 80 STR 
LOGICAL*! NEGATIVE , DONE 

READ (5, 1000) STR 
1000 FORMAT (A80) 

K - 0 
I - 1 

DONE - .FALSE. 

^J^ofo^*" 07 * °° NE) • AND - (K - LT - H0W >> 

DO WILE ((I ; LT. 80) .AND. (STR(I:I) .HE. '-•> .AND. 
J (STR(I:I) ,NE. ' + •) .AND. (STTtfT *T\ m t ■ % 

x ^ ((STR(I:I) .LT. . 0 -J .OR?* C if5?if ij> . fiff • . ^ r J j , 

END DO 

IF (I .LT. 80) THEN 

x ^J5I51l: I) - EQ - then 

NEGATIVE - .TRUE. 
I - I + 1 

EL 5£-?f J STR < I:I > -EQ. '*•) THEN 

NEGATIVE - .FALSE. 

1-1*1 
ELSE 

NEGATIVE - .FALSE. 
END IF 

r Vl°l^ f r (is ' s '^ s:>,: ^*^ SR(1!I) ,LE ,9,)) 

00"^\(I^W. I 80,^AKO. (STR (I!1) . CE . 'Q • ) •AND. 

i - r + i° * <1CHAR < STR < 1J i>> - «> 

END DO 

IF^(STR(I:I) . EQ. •••) THEN 

% ""™!ifi f xS> ?&^V STR(rsI) - GE - ,Q,) •*»• 

1-1*1°* f ICHAR ( STR (I-D) - 48) 
END DO 
M - 2 
J2 - J 

DO WHILE (M .LE. I-L) 
J2 - J2 / 10 
M - K ♦ 1 
END DO 
END IF 



WO 90/05920 



PCT/US89/04945 



66 



ELSE 

I * tt fW«*TtVI) THEN 
hJ^SW * J2 * (-1, 

HDMS(K) - J2 
BHD IF - ^* 
ELSE 

END DO 

RETDRR- 
END 

FUNCTION GET r 
Jin*GERM CEf t 

IOGICAL*! NEGATIVE 



J - 0 



% °° f(t 80) -AND. (STRU.-I, ,NE. ^ 

* _ ■ . (fSTRfi: Z ) .lt. {a 2 (t 52i-:"«i. .'and." 

BHD DO 1 CSTR(I:I) .ct. '9^)) 

X * X x ' 
ZXSE 

^NEGATIVE - .FAI5E.. 

I -X * l° * / ICH ^CSTR C X:I)) i } 48) 
END DO 
END IT 

XF (NEGATIVE) THEN 
END IF * 
GET I - j 

retOrn 

END 
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::::::::::::::::::.s^ 

! 

f COPYRIGHT - KAYO FOUNDATION FOR MEDICAL EDUCATION AND RESEARCH - 198 8 

SUBROUTINE MAKB_NEW_UDF ( FN , IER4 ) 

INCLUDE , :RDISK:UTIL:QSYM.F77.IN' 

CHARACTER* (•) FN 
INTEGER* 4 IER4 

INTEGER* 2 CPACK (0: ISYS^CLTH) 

INTEGER* 4 ACO, AC1 # AC2 

JNTEGERM CTIK4, CACP4 

EQUIVALENCE (CTIM4 , CPACK ( I SYS CTIM) ) 

EQUIVALENCE (CACP4 r CPACK (IS YS^CACT) ) 

ACO - BYTEADDR ( FN) 

IER4 - ISYS(ISYS_DELETE,AC0,AC1,AC2) 

^gSS-SR' : 5 S "-^*256 ♦ ISYS.FUDF 
CTIM4 - 
CACP4 - -1 

CPACX(ISY5 CDEH) - 0 

CPACK (ISYS COEL) - -1 

CPACK(ISYS CMIL) - -1 

CPACK ( IS YS^OfRS) - 0 

ACO - BYTEADDR ( FN ) 
AC2 - WORDADDR ( CPACK) 

IER4 - ISYS(ISYS_CREAT£,AC0,AC1,AC2) 

RETURN 
END 

rnilll!!!!!!!!*?^*^ 

»*or make new_udf. F77 source file 
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Claims 

We claim: 

1. An NMR system, the' combination comprising: 

means for generating a polarizing magnetic field; 

excitation means for generating an RF excitation 
magnetic field which produces transverse magnetization 
in nuclei subjected to the polarizing magnetic field; 

receiver means for sensing the NMR signal produced 
by the transverse magnetization and producing digitized 
in-phase (I) and. quadrature (Q) samples of the NMR 
signal; 

first gradient means for generating a first 
magnetic field gradient to impart a first phase 
component into the. NMR signal which is indicative of the 
location of the transversely magnetized nuclei along a 
first coordinate axis; 

second gradient means for generating, a second 
magnetic field, gradient to impart a second phase 
component into the NMR signal which is indicative of the 
location of the transversely magnetized nuclei, along a 
second coordinate axis; 
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20 pulse control means coupled to the excitation 

means, first and second gradient means, and receiver 
means/ said pulse control means being operable to 
conduct a scan in which a series of pulse sequences are 
conducted in which the second magnetic field gradient is 

25 stepped through a series of discrete values and a 
corresponding series of NMR signals are sensed and 
digitized to form an NMR data set; and 

processor means for storing the NMR data set and 
for reconstructing an image array for a display from the 

30 stored NMR data set by: 

(a) Fourier transforming the NMR data set along 
one of its dimensions to create hybrid-space data arrays 
I" and Q f ; 

(b) producing a correction data array using the 
35 data in the hybrid-space data arrays I* and Q f ; 

(c) applying the data in the correction data array 
to the NMR data set to reduce motion effects; and 

(d) Fourier transforming the corrected NMR image 
data set to produce the image array. 

2. The NMR system as recited in claim 1 in which 
the correction data array is produced by calculating the 
magnitude of the transformed sampled NMR signals in the 
hybrid-space data arrays I' and Q 1 to produce a modulus 
5 array, and correlating each row of the modulus array to 
produce a corresponding shift value for the correction 
data array. 
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3 V The NMR system as recited in claim 2 in which 
a rollback correction data array is produced from the 
corresponding shift values in the correction data array 
and . a rollover rate Kr, and the rollback correction data 
array values are applied to correct the NMR image data 
set in step . (c> . 

4. The- NMR system as recited in claim 1 in which 
the correction data array is produced by calculating the 
phase of the transformed sampled NMR signals in the 
hybrid-space data arrays I ' and Q* to produce a two- 
dimensional phase array, and determining the difference 
in phase between elements in a reference row of the 
phase array and elements in the same column of the phase 
array to produce the values for the correction data 
array. 

.5- The NMR system as recited in claim 2 in which 
a second correction data array. is produced by 
calculating the phase of the transformed sampled NMR 
signals in the hybrid-space data arrays I' and Q* to 
produce a two-dimensional phase array, and determining 
the phase difference between elements in a reference row - 
of the phase , array and elements in the same column of 
the phase array to produce the values for the second 
correction data array. 
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6. The NMR system as recited in claim 1 in which 
the NMR data set acquired during the scan includes NMR 
navigator data which has been sub jected to one of said 
two magnetic field gradients and NMR image data which 

5 has been subjected to both of said two magnetic field 
gradients, and steps (a) and (b) are performed with the 
NMR navigator data and step (c) is performed on the NMR 
image data. 

7. In an NMR system which performs a scan to 
acquire an NMR data set from which an image array is 
reconstructed, the improvement comprising: 

transforming an NMR data set produced by the NMR 
5 system to create a hybrid- space data array; 

producing a correction data array using the data in 
the hybrid-space data array; 

applying the data in the correction data array to 
an NMR data set produced by the NMR system to reduce 
10 artifacts caused by motion in the acquired NMR data set; 
and 

producing an image array from the corrected NMR 
data set. 



8. The improvement as recited in claim 7 in which 
the hybrid- space data array is created by: 

performing a Fourier transformation on the NMR data 
set along one of its dimensions. 
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9. The improvement as recited in claim 8 in which 
the correction data array is produced by calculating the 
magnitude of each element of the hybrid-space data array 
to produce a corresponding modulus array, and 
correlating each row of the modulus array to produce a ' 
corresponding shift value for the correction data array. 

.10. The improvement as recited in claim 9 in which 
a rollback correction data array is produced from the 
corresponding shift values in the correction data array 
and a rollover rate K R , and the values in both the 
correction data array and the rollback correction data 
array are applied to correct said acquired NMR data set 
used to produce the image array. 

11. The improvement as recited in claim 8 in which 
the correction data array is produced by calculating the 
Phase of each element of the hybrid-space data array to 
produce the corresponding elements of a two-dimensional 
phase data array, and determining the difference in 
phase between elements in a reference row of the phase 
data array and elements in the same column of the phase 
data array. 
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12. The improvement as recited in claim 9 in which 
a second correction data array is produced by: 

calculating the phase of each element of the 
hybrid-space data array to produce the corresponding 
elements of a two-dimensional phase data array having 
rows and columns; and 

calculating the elements of the second correction 
data array by determining the difference in phase 
between elements in a reference row of the phase data 
array and elements in the same column of the phase data 
array; 

wherein the data in both the correction data array 
and the second correction data array is applied to the 
acquired NMR data set to reduce the effects of motion, 

13. The improvement as recited in claim 7 in which 
the NMR data set used to create the hybrid-space data 
array is acquired from a navigator NMR signal which is 
produced during each pulse sequence of a scan but which 
is not subjected to a changing magnetic field gradient 
during the scan, and the NMR data set employed to 
produce the image array is acquired from an NMR signal 
which is produced during each pulse sequence of the scan 
and which is subject to a phase encoding magnetic field 
gradient which changes during the scan. 
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14. - The improvement as recited in claim 13 in 
which each navigator NMR signal is produced in the same 
pulse sequence with its corresponding phase encoded NMR 
signal. 

15. An NMR system, the combination comprising: 
means for generating a polarizing magnetic field; 
excitation means for generating an RF excitation 

magnetic field which produces transverse magnetization 
in nuclei;, subjected to the polarizing magnetic fields- 
receiver means for sensing the NMR signal produced 
by the transverse magnetization and producing digitized 
in-phase (I) and quadrature (Q) samples of the NMR 
signal; 

first gradient means for generating a first 
magnetic field gradient to impart a first phase 
component into the NMR signal which is indicative of the 
location of the transversely magnetized nuclei along a 
first coordinate axis; 

second gradient means for generating a second 
magnetic field gradient to impart a second phase 
component into the NMR signal which is indicative of the 
location of the transversely magnetized nuclei along a 
second coordinate axis; 
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20 pulse control means coupled to the excitation 

means, first and second gradient means, and receiver 
means, said pulse control means being operable to 
conduct a scan in which a series of pulse sequences are 
conducted in which the second magnetic field gradient is 

25 stepped through a series of discrete values and a 
corresponding series of NMR signals are sensed and 
digitized to form an NMR data set; and 

processor means for storing the NMR data set and 
for reconstructing an image array for a display from the 

30 stored NMR data set by: 

(a) producing a phase array from the NMR data set 
which indicates the phase of the digitized NMR signals; 

(b) producing correction data using the data in 
the phase array; 

35 (c) applying the correction data to the NMR data 

set to reduce motion effects; and 

(d) Fourier transforming the corrected NMR image 
data set to produce the image array. 
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16.. An NMR system, the combination comprising: 

means for generating a polarizing magnetic field; 

excitation means for generating an RF excitation 
magnetic field which produces' transverse magnetization 
in nuclei subjected to the polarizing magnetic field; 

receiver means for sensing the' NMR signal produced 
by the transverse magnetization and producing digitized 
in-phase (I) and quadrature (Q) samples of the NMR 
signal; 

first gradient means for generating a first, 
magnetic field gradient to impart a first phase 
component into the NMR signal which is indicative of the 
location of the transversely magnetized nuclei along a 
first coordinate axis; 

second gradient means for generating a second 
magnetic field gradient to. impart a second phase 
component into the NMR signal which is indicative of the 
location of the transversely magnetized nuclei along a 
second coordinate axis; 

pulse control means coupled to the excitation 
means, first and second gradient means, and receiver 
means, said pulse control, means being operable to 
conduct a scan in which a series of pulse sequences are 
conducted in which the. second magnetic field gradient is 
stepped through a series of discrete values and a 
corresponding series of NMR signals are sensed and 
digitized to form an NMR data set; and 
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processor means for storing the NMR data set and 
for reconstructing an image array for a display from the 
30 stored NMR data set by: 

(a) producing a phase array from the NMR data set 
which indicates the phase of the digitized NMR signals; 

(b) producing an altered phase array by 
subtracting phase components produced by stationary 

35 spins from each value in the phase array; 

(c) producing correction data using the data in 
the altered phase array; 

(d) applying the correction data to the NMR data 
set to reduce motion effects; and 

40 (e) transforming the corrected NMR image data set 

to produce the image array. 

17 . The NMR system as recited in claim 16 in which 
the phase components produced by stationary spin are 
calculated by determining the center of a circle which 
best fits the values in each column of the phase array. 
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